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COMPUTER PROGRAM FOR POST-FLIGHT EVALUATION OF THE CONTROL 
SURFACE RESPONSE FOR AN ATTITUDE CONTROLLED MISSILE 

SUMMARY 


A FORTRAN IV coded computer program is presented for post-flight analysis 
of a missile's control surface response. It Includes preprocessing of 
digitized telemetry data for time lags, biases, non-l-inear calibration changes 
and filtering. Measurements include autopilot attitude rate and displacement 
gyro output and four control surface deflections. Simple first order lags are 
assumed for the pitch, yaw and roll axes of control. Each actuator is also 
assumed to be represented by a first order lag. Mixing of pitch, yaw and roll 
commands to four control surfaces is assumed. A pseudo-inverse technique is 
used to obtain the pitch, yaw and roll components from the four measured 
deflections. 

This program has been used for over 10 years on the NASA/SCOUT launch 
vehicle for post-flight analysis and was helpful in detecting incipient 
actuator stall due to excessive hinge moments. 

The program is currently set up for a CDC CYBER 175 computer system. It 
requires 34K words of memoiy and contains 675 cards. A sample problem 
presented herein including the optional plotting requires eleven (ll) seconds 
of central processor time. 
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1.0 INTRODUCTION 


Post-flight analysis of a missile autopilot and control system should 
include comparison of measured parameters for consistency with preflight 
mathematical models. For a proportional control system telemetered data may 
include the autopilot gyro outputs and the measured control surface 
deflections. Preflight gains applied to the gyro outputs and the model should 
yield the control surface deflections. Differences between the reconstructed 
values and measured values may be indicative of anomalies or incipient 
failures. This report presents a computer program used for post-flight 
analysis of the NASA/SCOUT Launch Vehicle first stage proportional control 
system. It has been used to identify anomalous behavior. 

The missile autopilot is assumed to contain three axes of information each 
of which includes a displacement and rate term such as, pitch, yaw and roll 
axes having angular displacement and rate included in the control law. A 
block diagram of the system is presented in Figure 1. 

In post-flight data reduction and analysis there are sometimes small 
deviations in time for each parameter, differences between assumed linear 
calibration for data reduction and actual non-linear calibrations, 
cross-coupling between on-board telemetry channels, and, higher frequency data 
and noise. All of these can have a significant effect on the post-flight 
reconstruction process. Therefore, a large part of the methodology herein 
involves shifting, adjusting and smoothing of the digitized reduced telemetry 
data. 


The assumptions, methodology, program description and running instructions 
are presented in the following sections. 
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2.0 METHODOLOGY 


This section conteiins the methodologj and equations which are used to 
adjust the telemetry data and to reconstruct the control surface response from 
the measured data. Telemetry data for the autopilot pitch, yaw, and roll 
channel attitude rates and displacements and four control surface deflections 
are required. Estimates or preflight measurements of autopilot gains are also 
required. A proportional control system represented hy the block diagram of 
Figure 1 is assumed. The assumptions and equations are presented in the 
following paragraphs. 


2.1 Assumptions 


Major assumptions and approximations are: 

. autopilot and control system as presented in Figure 1, 

. control gains are constant in time, 

. pitch, yaw, and roll channels frequency response can be represented as 
a first order lag (single break frequency), 

. actuator response for each of the four control surfaces can be 
represented by a first order lag, 

. mixing of the pitch, yaw, and roll error signals for commands to the 
four control surfaces is represented by a constant matrix 
transformation, 

. frequency response of the telemetry data is greater than the autopilot 
model 

. phase shifts in the telemetry data can be represented by an equivalent 
time shift for each of the parameters, 

. telemetry system cross-coupling is limited to the attitude displacement 
gyro output and is proportional to vehicle rate about that axis 

. non-linearities in calibration of the control surfaces are not time 
dependent 
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2.2 Equations 


2.2.1 Control Surface Response 


The control surface response is represented by the block diagram of Figure 
1. The pitch, yaw, and roll channel error signals are. 


(2.1) 



(2.2) 



(2.3) 


+ I^R^P 


+ 


These error signals are processed through a linear transformation to mix 
the signals to each of four control surfaces. This transformation is. 


(2.4) 


»c- W 


where. 


is the four element vector 


(2.5) 


Sc = «C2 


[A] is the 4 by 3 transformation matrix specified by the autopilot design 

(2.6) [a] - rail ®12 


a2i a22 a23 

®31 ®32 *33 

*41 *42 *43_ 

is the three element error vector 


According to Figure 1 each actuator responds to commands as a simple first 
order lag, i.e. , , 

«•( **act ) 


8 = Sr 


s + (t) 
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Up to this point the system is linear. It is assummed that each actuator 
responds with the same break frequency. Therefore, the transfer function of 
equation (2.7) can be moved and included in the transfer function of equations 
(2.1), (2.2), and (2.3)* This is done in the program to reduce computer 
time. The commanded pitch, yaw, and roll surface deflections become. 



This will be followed through the derivation for the pitch channel; yaw 
and roll channel equations are similar. Prior to the mixing of channel error 
signals via the [aJ transformation, the equivalent pitch component of 
commanded surface deflections is. 



This can be put into the state variable form. 


(2.10) X = [f] X + b u 


where the iiqjut. 



(2.14) b = 


and the output is, 
(2.15) X2 


(two state vector) 

(system matrix) 

(input coefficient vector) 


In the program these equations are solved by a linear system time response 
subroutine (TRESP) which uses a fourth-order RUNGE-KUTTA integration method. 
After the pitch, yaw, and roll commanded deflections are computed they are 
transformed by the matrix [a] (Equation 2.6) to provide the calculated or 
reconstructed deflections of each of the four control surfaces. 


2.2.2 Pseudo-Inverse 


Commanded individual control surfaces represented by equation 2.4 involves 
a linear transformation [a] of a 'three-vector' into a 'four-vector'. The 
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pitch, yaw, and roll components of the four control surfaces can be calculated 
from the measured individual control surface deflections* This is easily done 
using the pseudo-inverse transformation of [a]. A matrix inverse does not 
exist for non-square matrices. The pseudo-inverse provides a reverse 
transfomation in a least-squares sense. The pseudo-inverse is, 

(2.16) [A]# - ([a]T[a])-1 [aP 

for the overdetermined case where [a] has more rows than columns. 

Reference ’’Applied Optimal Estimation" edited by Arthur Gelb, published by 
the M.I.T. Press, Cambridge, Massachusetts, 1974. 

Actual pitch, yaw, and roll control surface deflection components can be 
computed, -| 

(2.17) 

2.2^3 Telemetered Data Adjustments 


’act^ 
^act,^ 
L^act^J 




"act]^ 

^act2 

^act3 

L^act4 


Reduced telemetiy data for the pitch, yaw, and roll displacements and 
rates and the four control surface deflections may require further adjustments 
and filtering. These adjustments include biases, time shifts, 
non-linearities, cross-coupling, and additional filtering. 

BIASES 


Each parameter is assumed to have a simple bias error which can be 
estimated from a quiet period of flight such as prior to vehicle liftoff. 
These are entered in the input to the routine and added to the reduced 
telemetiy data, i.e., 

(2.18) yj^(t) = y^(t) +^y 

where, ym(t) is the reduced measured telemetry parameter 
y5(t) is the adjusted parameter 
Ay is the bias shift required 

TIME SHIFTS 


Each measured parameter is assumed to have a different time delay due to 
the telemetry system and ground station playbacks. In order to be consistent 
each parameter must be shifted to a common time base. This is done by a table 
lookup versus time such that the parameter (y) is shifted back in time by its 
peculiar lag, i.e., 

(2.19) y"ni(t) = y’m(t+At) 

where, y"ni('*') ia "the measured parameter adjusted for time and bias 

shift s 

is the required time shift for parameter 'y' 


- 5 - 


At 



NON-LINEARITIES 


The four control surface deflections may have a parabolic calibration 
curve (this is the case for a measurement using a linear potentiometer 
attached to a bellcrank). In data reduction the calibration applied is 
usually a straight line or a series of straight lines (see the sketch below). 




This can lead to subtle errors in reconstructing the control surface 
deflections. Therefore, an adjustment is made to the data via a table lookup 
in the computer program. This adjustment is, 

(2.20) Ay = Ay-biag + Ay (5) 

where, Ayv^og is average bias value 

AyCo) is a function of the deflection based on the sketches 
above 

The bias shown in Equation (2.20) is inserted into Equation (2.18). 

CROSS- COUPLING 


Peculiarities in an attitude gyro telemetry pickoff may induce a voltage 
proportional to the vehicle rate about the input axis (this is peculiar to the 
SCOUT vehicle using miniature integrating rate gyros to measure attitude 
displacements). The true attitude displacement is therefore. 


( 2 . 21 ) 


00 

</fe 


, m 

K- " 
% - p 


where the coefficients a^, and a^ are determined by preflight test of 
the telemetry and autopilot subsystems. 


ADDITIONAL FILTERING 


Since the telemetry data may have noise or data at frequencies well above 
the bandwidth of interest, additional filtering capability is afforded by the 
computer program. A third order Butterworth filter is supplied. The cutoff 
frequency is an input variable. Each parameter is filtered with the same 
algorithm and cutoff frequency. The transfer function is, 

(2.22) T(s) = 

1 + 2 r s + 2 
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where, 

(2.25) r= 0.707 /cOj,Q 
aod, 



state is the filter output 

These equations are solved in the computer program using the time response 
subroutine ( TRESP) . 


Application of this filter adds a low frequency time delay to the data. 
This delay is, 

(2.25) t = 1.418 A>qo (seconds) 

This time increment is added to the individual telemetry parameter time 
delays before shifting the data with Equation (2.19). 
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5.0 PROGRAM DESCRIPTION 


3.1 General 


This computer program is programmed in FORTRAN IV for a CDC CYBER 175 
system. The coding is in the most part compatible with ANSI standards. 
Non-ANSI statements include the PROGRAM card and the use of EOF (end-of-file) 
test for transfer from reading input. Another area of limited portability is 
the use of ten letter words for labeling information. Consult your computer 
department for changes in these areas. 

The computer routine is arranged to operate with standard card input and 
line printer output. Optional plotting is based on standard CALCOMP plotters 
and software. The CALTEKA Tektronix terminal plot can also be used without 
program modification. 

A main routine (FINRES) and twelve (12) subroutines require approximately 
34K words of computer memoiy. Input and output is stored in arrays to 
facilitate time adjustments, filtering, and a well formatted printed and 
plotted output. 

Program flow and user instructions are presented in the following 
paragraphs. Input and output of a sample problem is illustrated along with 
detailed descriptions. 


3.2 Program Flow 


Program flow is straightforward in eight basic parts. 

. input data 

. telemetry data adjust for calibration and filtering 
. reconstruction of pitch, yaw and roll commands 

. comparison of individual control surfaces with reconstructed commands 
. output of individual surface data 

. optional reconstruction of pitch, yaw and roll components of control 
surface response and output 
. plotted output 

The interdependence of the main routine and the subroutines is presented 
in Figure 2. A flow chart of the main routine (FINRES) is presented in Figure 
3. A complete listing of the FORTRAN program and subroutines other than the 
standard CALCOMP library subroutines are presented in Appendix A. 

Descriptions of the twelve subroutines are presented in the following 
paragraphs. 


- 8 - 



3*3 Subroutine Description 


Twelve subroutines are used to support the FINRES main program; CURVE, 
DASH, ERSIG, FILFIL, FIN, PSEUDO, RUNGE, SIMEQ, TBLN, TRESP, XMULT, and YDOT. 
A brief description of each is presented below. 

CURVE 


This subroutine sets up the calcomp plot for one frame of a single control 
surface deflection comparison. It includes the graph paper description 
(CAL22) which has a 10 by 16 grid size with 20 divisions per inch. Paper 
performation size is 11 by 17. All data to be plotted by this routine enters 
in the argument list. The actual curve plotting is made through calls to the 
DASH subroutine. The call statement for CURVE is, 

CALL CURVE (T, CALC, ACT, NP, NTIT, NM, XS, YS, DS) 

where , 

T input array name of time abscissas 

CALC input array name of computed commanded surface deflections 

(reconstructed) 

ACT input array name of measured control surface deflection 

NP is the number of time points in T, CALC, and ACT arrays to 

be plotted 

NTIT input array of eight (8) ten-letter words contained 80 

character title 

NM input ten letter word variable to identify frame (this is 

output on second line of title) 

XS abscissa (time) scale factor (units/inch) 

YS ordinate deflection scale factor (units/inch) 

DS ordinate scale factor for difference (CALC-ACT) deflection 

(units/inch) 

Care must be taken in selecting scale factors so that plotted data falls 
on grid. Limiting of plotted data is automatically invoked in CURVE through 
the call statements to DASH. 


DASH 

This subroutine plots a curve on a CALCOMP plotter for a set of 
ordinates and abscissas. The style and type of line drawn is selected by the 
user. Note that the CALCOMP plot is specified in inches; plotting on metric 
paper requires appropriate scaling change before entering this subroutine. 

The call statement is, 

CALL DASH (X, Y, NP, Zl, Z2, SPACE, XSCALE, YSCALE, LSYMB, XLIM, YLIM) 
where, 

X - input array of abscissa values 

Y - input array of ordinate values 

NP - number of points in X and Y to be plotted 
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Z1 - for dashed-dot lines this is length, of long line measured 

in inches’ (see sketch helow) 

Z2 - for dashed-dot lines this is length of short line measured 

in inches (see sketch) 

SPACE - for dashed style lines this is the length of the space between 
lines measured in inches. 

SPACE = 0 gives a solid line plot 

SPACE =■ negative gives special CALCOMP symbols at each point 
XSCALE - abscissa plot scale factor (units per inch) 

YSCALE - ordinate plot scale factor (units per inch) 

LSYMB - special CALCOMP symbol code number used if SPACE is negative 
(see code below) 

- (+) I5YMB gives straight solid lines between symbol points 

- (-) LSIMB gives only symbols at each point without lines 
XLIM - plot limiting of the abscissa (inches) points out of range, 

range will appear at this limit 
YLIM - plot range of ordinate (inches) 

For ease in use, the following styles are typically possible. 



TYPE 

Z1 

Z2 

SPACE 

LSYMB 

Solid 

— 

— 

0. 

0. 

Dashed 

0.25 

0.25 

0.10 

0. 

Dashed 

0.07 

0.07 

0.07 

0. 

Dashed 

Dot 

0.5 

0.03 

0.07 

0. 

Sjrmbols 

— 

— 

-0.1 

+2 

Symbols 
(no line) 

— 

— 

-0.1 

-2 


ERSIG 


This subroutine computes the single axis error signal commanded deflection 
(pitch, yaw, or roll). It includes adjustment for rate crosscoupling into the 
displacement telemetry data, the channel break frequency and actuator break 
frequency time response. The call statement is, 

CALL ERSIG (FILTER, B, T, Q, TH, NP, NINT, CTH, CTHD, DT, DTTM, DTTMR, AD, 
V, WACT) 

where, 

FILTER input third order coefficient array for the Butterworth 

filter states (see Paragraph 2.2.5) 

B telemetry filtering input coefficient vector (Paragraph 

2.2.3) 

T input array of time points for telemetry data containing NP 

points 

Q input array of NP points of telemetered rate data 

corresponding to (T) times. It is changed to output the 
filtered error signal. 
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TH input array of NP points of telemetered displacement data 

corresponding to (T) times 

NP number of time points in arrays T, Q, and TH 

HINT niuaber of integration steps per value of time in (T) array 

to be used in filtering and response histories (e.g., if 
HINT = 2) the integration step size is one-half of the time 
between points in (t) array 
CTH attitude displacment gain Kp 

CTHD rate gain Kp 

DT time delay of third order filter (FILTER) and (B) to be used 

in time corrections 

DTTM time delay of attitude error telemetry data 

DTTHR time delay of rate telemetry data 

AD telemetry cross-coupling coefficient of attitude 

displacement due to rate 

W break frequency of the control channel 

WACT actuator break frequency 

Note that the array (Q) is destroyed by the subroutine and used to return 
the computed filtered error signal. 


FILFIL 


This subroutine fills the third order Butterworth filter coefficient 
arrays and computes the low frequency time lag of this filter as presented in 
Paragraph 2.2.3* The call statement is, 

CALL FILFIL (WCO, A, B, DT) 
where, 

WCO input cutoff frequency of the Butterworth filter (radians 

per second) 

A output coefficient matrix of the filter (3 by 3) 

B output coefficient vector for filter state inputs 

DT output effective low frequency time lag associated with the 

filter (seconds) 

FIN 

This subroutine filters a telemetered control surface deflection time 
history and adjusts it for time shifts. The call statement is, 

CALL FIN (FILTER, B, T, D, NP, DTFIL, TC, NINT) 
where, 

FILTER input (3 by 3) matrix of third order filter coefficient 

matrix 

B input (3 by l) coefficient input vector for third order 

filter 

T input array of NP times for deflection array D 

D input array of NP control surface deflections (also the 

output filtered time adjusted deflection) 

NP number of time points in T, and D arrays 

DTFIL input low frequency time lag of filter in seconds 

TC input time lag of telemetered control surface deflection for 

use in adjustment 
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HINT 


number of integration steps to be used between time points 
in (T) array 


PSEUDO 


This subroutine computes the pseudo - inverse of the coefficient matrix as 
described in Paragraph 2.2.2. The call statement is, 

CALL PSEUDO (B, A, N, M, NER) 

where, 

B output M by N pseudo - inverse matrix of A 

A input matrix having dimensions N by M 

N is number of rows of (A) and number of coliimns of (B) 

M is number of columns of (A) and number of rows of (B) 

NER 13 an error indicator 

NER = 1 nomal execution 

NER = 0 abnormal condition because of a submatrix 

singularity (pseudo inverse cannot be computed) 

RUNGE 


This subroutine aids in the RUNGE-KUTTA integration of the time response 
(TRESP) subroutine. The call statement is, 

CALL RUNGE (N, FN, H, X, Y, L, l) 

where, 

N input system order 

FN first derivatives of the state 

H integration step size 

X time variable 

Y state vector 

L control integer 

L = 1 indicates incomplete integration process 
L = 2 indicates completed integration step 
I number of times that RUNGE has been entered on the current 

integration step (when 1=5 the final answer is computed) 

SIMEQ 

This subroutine is used to compute the inverse of a square matrix by 
diagonalization. With modification it can be used to solve a set of 
simultaneous linear equations. The call statement is, 

CALL SIMEQ (A, KC, AINV, lERR) 

where, 

A input KC by KC matrix to be inverted 

KC input order of the matrix 

AINV output inverted matrix if computed 

lERR output error variable 

lERR = 1 normal computation 

lERR - 0 abnormal (A matrix is singular) 
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TBLN 


This is a single table lookup subroutine using linear interpolation 
between points. This subroutine requires separate consistent arrays of 
abscissas and ordinates. The abscissas must be in ascending order. The call 
statement is, 

CALL TBLN (Y, X, T, A, NT, M) 
where, 

y output ordinate to be found 

X input abscissa value 

T array of abscissas 

A array of ordinates 

NT number of values in (T) and (A) arrays 

M input index 1 to NT to begin the table search. After 

locating the ordinate the nearest location is returned for 
future use 


TEESP 


This subroutine used in conjunction with subroutines RUNGE and YDOT solve 
a time response for a set of up to three first order linear differential 
equations by a fourth order RUNGE-KUTTA integration procedure. The call 
statement is, 

CALL TRESP (A, B, T, Y, Z, NP, N, K, NINT) 
where, 

A input coefficient matrix 

B input coefficient vector for state equations 

T input array of (NP) forcing function time values 

Y input array of NP forcing function values corresponding to 

times in (T) array 

Z output array of 'Kth' state variable values corresponding to 

time in (T) array 

NP input niunber of values in T, Y, and Z array 

N input order of the system 

K input designation of state variable corresponding to the 

desired output 

NINT input number of integration steps between time points in (T) 

array 


XMULT 


This subroutine multiplies two matrices. The call statement is, 
CALL XMULT (A, B, C, N) 
where, 

A input premultiplier matrix 

B input postmultiplier matrix 
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C output matrix product of (A.B) having order N by N 

N is desired order of the output matrix (if input matrices are 

non-square they are filled with zeroes where rows or columns 
are needed) 


YBOT 

This subroutine computes the first derivative of the system state vector 
used by TRESP to integrate a set of first order linear differential 
equations. The call statement is, 

CALL YDOT (A, Y, XDOT, B, U, N) 
where, 

A input system coefficient matrix of order N 

Y input state vector 

XDOT input derivative of the state vector and returned updated 

derivative of the state vector 

B input coefficient vector of the forcing function 

U input value of the forcing function 

N input order of the system 


3.4 Input Data Description 


Input data descriptions are presented in the following subparagraphs. A 
sample problem input data listing is presented in Figure 4 for reference. 
Input data can be separated into the following categories; 

1) title card and control options 

2) telemetry adjustements 

3) control sytem constants 

4) control surface deflection adjustment tables for non-linearities 

5) plot scale factors 

6) time histories of measured telemetry variables 


3.4.1 Title and Control Constants 


The first card of input contains 80 columns of arbitrary title 
information. This information is printed at the head of each page of output 
and at the top of each frame of plotted data. It is input using an array 
(NTIT) having eight words containing ten characters each. 

The second card contains four (4) integer constants right justified 
without a decimal point. It is input with a format of (415). The description 
is. 
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FORTRAN 

NAME 


COLUMN 


DESCRIPTION 


IOC 


lOP 


NPRT 


NINT 


5 control option 

IOC =0 do not compute pitch yaw and 

roll deflection components using 
pseudo-inverse 

IOC =• 1 compute pitch, yaw, and roll 
deflection components using 
pseudo-inverse 

10 plot option 

lOP =■ 0 no plots 

lOP = 1 plot time histories of control 
surface deflections 

11-15 printed output control integer NPRT + 1 

time intervals are skipped in output 
(e.g., NPRT = 2 keeps every other point 
from being printed) 

16-20 number of integration steps for each 

increment of time between the equally 
spaced input data 


3.4*2 Telemetry Adjustments 


This group contains twenty-four (24) constants used for adjustment of the 
telemetry data. These are input eight numbers per card in fields of ten 
columns (format 8E10.3). The descriptions are. 


FORTRAN 

NAME 

CARD 

COLUMNS 

DESCRIPTION 

DL(1) 

3 

1-10 

bias adjustment to be added to pitch 
rate 

DL(2) 

3 

11-20 

bias adjustment to be added to yaw rate 

DL(3) 

3 

21-30 

bias adjustment to be added to roll rate 

DL(4) 

3 

31-40 

bias adjustment to be added to pitch 
displacement 

DL(5) 

3 

41-50 

bias adjustment to be added to yaw 
displacement 

DL(6) 

3 

51-60 

bias adjustment ot be added to roll 
displacement 

DL(7) 

3 

61-70 

bias adjdustment to be added to fin 1 
control surface 
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FORTHAN 


NAME 

CARD 

COLUMNS 

DESCRIPTION 

DL(8) 

3 

71-80 

bias adjustment to be added to fin 2 
control surface 

DL(9) 

4 

1-10 

bias adjustment to be added to fin 3 
control surface 

DL(IO) 

4 

11-20 

bias adjustment to be added to fin 4 
control surface 

TCQ 

4 

21-30 

time lag in pitch rate data 

TCR 

4 

31-40 

time lag in yaw rate data 

TCP 

4 

41-50 

time lag in roll rate data 

TCTH 

4 

51-60 

time lag in pitch displacement data 

TCPS 

4 

61-70 

time lag in yaw displacement data 

TCPH 

4 

71-80 

time lag in roll displacement data 

TCI 

5 

1-10 

time lag in fin 1 control surface data 

TC2 

5 

11-20 

time lag in fin 2 control surface data 

TC3 

5 

21-30 

time lag in fin 3 control surface data 

TC4 

5 

31-40 

time lag in fin 4 control surface data 

AKTH 

5 

41-50 

pitch displacement to rate telemetry 
cross coupling coefficient 

AKPS 

5 

51-60 

yaw displacement to rate telemetry 
cross coupling coefficient 

AKPH 

5 

61-70 

roll displacement to rate telemetry 
cross coupling coefficient 

WCO 

5 

71-80 

cutoff frequency of third order 
Butterworth filter to be used on all 
input time histories (hertz) 


3 *4. 5 Control System Constants 


This group of input data includes the autopilot control system gains, 
break frequencies and error signal mixing matrix for the individual control 
surface commands. Gains and break frequencies are input in fields of ten with 
format (8E10.3). They are, 
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FORTRAN 

NAME 

CARD 

COLUMN 

SYMBOL 

UNITS 

DESCRIPTION 

CTH 

6 

1-10 

Kd^ 

— 

pitch attitude 
displacement gain 

CTHD 

6 

11-20 

Kr0 

sec 

pitch rate gain 

CPS 

6 

21-50 

^D,/r 

— 

yaw attitude displacement 
gain 

CPSD 

6 

31-40 


sec 

yaw rate gain 

CPH 

6 

41-50 

% 

— 

roll attitude 
displacement gain 

CPHD 

6 

51-60 

Kr0 

sec 

roll rate gain 

VIQ 

6 

61-70 


rad/sec 

pitch channel break 
frequency 

¥1R 

6 

71-80 

W,/r 

rad/sec 

yaw channel break 
frequency 

¥1P 

7 

1-10 


rad/sec 

roll channel break 
frequency 

¥ACT 

7 

11-20 

act 

rad /sec 

control surface actuator 
break frequency 


The next four cards contain the transfer matrix describing the 
relationship of the four control surfaces to the pitch, yaw and roll error 
signals. These are input with one card for each surface with three constants 
relating the components of pitch, yaw and roll error signal to each surface- 
input is with format (4E10.3). 

FORTRAN 


NAME 

CARD 

COLUMNS 



DESCRIPTION 

A(l,l) 

8 

1-10 

fraction 

of 

pitch error for surface no. 1 

A(l,2) 

8 

11-20 

fraction 

of 

yaw error for surface no. 1 

A(l,2) 

8 

21-30 

fraction 

of 

roll error for surface no. 1 

A(2,l) 

9 

1-10 

fraction 

of 

pitch error for surface no. 2 

A(2,2) 

9 

11-20 

fraction 

of 

yaw error for surface no. 2 

A(2,3) 

9 

21-30 

fraction 

of 

roll error for surface no. 2 

A(3,l) 

10 

1-10 

fraction 

of 

pitch error for surface no. 3 

A(3,2) 

10 

11-20 

fraction 

of 

yaw error for surface no. 3 

A(3,3) 

10 

21-30 

fraction 

of 

roll error for surface no. 3 

A(4,l) 

11 

1-10 

fraction 

of 

pitch error for surface no. 4 

A(4,2) 

11 

11-20 

fraction 

of 

yaw error for surface no. 4 

A(4,3) 

11 

21-30 

fraction 

of 

roll error for surface no. 4 
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3.4.4 Control Surface Non-Linear Calibration Adjustment Tables 


This group of input includes a table of adjustments for each control 
surface which allows for non-linearities in calibration not included xn the 
data reduction process. These are read with format (I5/»8E10.3) . The first 
card of each table contains the number of abscissa-ordinate pairs. The 
abscissa (control surface deflection) must be in ascending order. Four 
tables, one for each fin control surface includes, 

FORTRAN 

NAME COLUMN DESCRIPTION 


NTl 

1-5 

number of pairs of abscissas and ordinates 
in surface no. 1 adjustment table 

CDl(l), EDI(I), 

1-10, 11-20, etc 

CDl contains abscissa value of surface no. 
1 deflection 

EDI contains ordinate value of surface no. 
1 to be added to input telemetry data 

NTl 

CD2(I) ,ED2(I), 

1-5 

1-10, 11-20, 

(similar to above description for control 
surface no. 2) 

NT3 

CD3(I),BD3(I), 

1-5 

1-10, 11-20, 

(similar to above description for control 
surface no. 3) 

NT4 

CD4(I),ED4(I), 

1-5 

1-10,11-20 

(similar to above description for control 
surface no. 4) 


5 . 4.5 

Plot Constants 


This includes a single card of scale factors for the CALCOMP type plots. 
This card is input only when (lOP = l) on the second card of input (Paragraph 
3 . 4 . 1 ). Input uses format ( 8 EIO. 3 ). The descriptions are. 


FORTRAN 


NAME 

COLUMNS 

UNITS 

XSCALE 

1-10 

sec/inch 

YSCALE 

11-20 

deg/inch 

DSCALE 

21-30 

deg/inch 


DESCRIPTION 

time abscissa scale factor per inch of 
paper 

control surface ordinate scale factor 
per inch of plot paper 

control surface deflection difference 
scale factor per inch of paper 
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3.4.6 Telemetered Data Time Histories 


This group of data is entered with format (1IF7.3); each card contains all 
of the telemetered parameters for a given time point. The time points must be 
equally spaced. The analysis uses all the input data time points and stops 
reading input when an *End-of-File' card is read in the input data stream. 

The description of this data follows (refer to the sample problem input of 
Figure 4). 

FORTRAN 


NAME 

COLUMN 

SYMBOL 

UNITS 

DESCRIPTION 

T(J) 

1-7 

t 

seconds 

time 

Q(J) 

8-14 

P 

deg/sec 

telemetered pitch rate 

R(J) 

15-21 

q 

deg/sec 

telemetered yaw rate 

P(J) 

22-28 

r 

deg/sec 

telemetered roll rate 

TH(J) 

29-35 


degrees 

telemetered pitch 
displacement error 

PS(J) 

36-42 


degrees 

telemetered yaw 
displacement error 

PH(J) 

43-49 

0e 

degrees 

telemetered roll 
displacement error 

D1(J) 

50-56 

Si 

degrees 

telemetered fin 1 control 
surface deflection 

D2(J) 

57-63 

S2 

degrees 

telemetered fin 2 control 
surface deflection 

D3(J) 

64-70 

S3 

degrees 

telemetered fin 3 control 
surface deflection 

D4(J) 

71-77 

S4 

degrees 

telemetered fin 4 control 


surface deflection 


3.5 Output Data Description 


Output includes printed data and optional CALCOMP type plots (if lOP = 

1). A detailed description of the output is presented in the following 
paragraphs with a sample problem for reference. The printed output includes 
three parts depending on option. 

1) individual control surface response 

2) pitch yaw and roll control surface response (if IOC = l) 

3) CALCOMP plots (if lOP = l) 
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3.5.1 Individual Control Surface Response 


Printed output for a sample problem is presented in Figure 5. The first 
part of the output includes time histories of these parameters for each of the 
four control surfaces. These parameters are. 


OUTPUT 

LABEL 

SYMBOLS 

UNITS 

COMMAND 

Sc 

deg 

ACTUAL 

Sact 

deg 

DELTA 

^c-^act 

deg 


DESCRIPTION 

reconstructed control surface 
deflection based on the telemetered 
autopilot data 

filtered smoothed and adjusted 
telemetry data for the control surface 
deflection 

difference between the reconstructed 
and actual deflection data 


3 . 5.2 Pitch, Yaw, and Roll Component Response 


If the option (IOC = l) is chosen the pseudo-inverse of the control 
surface mixing gain matrix is computed. The pitch, yaw, and roll average 
components based on the four measured control surfaces are then computed. 
These are presented in Figure 6 for the sample problem. The parameters for 


the pitch, yaw, and roll channels are. 


OUTPUT 

LABEL 

SYMBOLS 

UNITS 

COMMAND 

Sc 

deg 

ACTUAL 

Sact 

deg 

DELTA 

Sc" ^act 

deg 


DESCRIPTION 

reconstructed pitch, yaw, and roll 
control error signal baaed on the 
telemetered autopilot data 

pitch, yaw, and roll component of the 
control surface deflection based on the 
telemetered control surface deflection 
and the pseudo-inverse gain matrix 

difference between the above values 


3 . 5.3 CALCOMP Plots 


If (lOP = 1 ) CALCOMP type plots are generated corresponding to the data 
printed. However, all computed points are plotted, whereas printed output can 
be suppressed by (NPRT) on the second input card. Sample problem plots are 
presented in Figure 7 and correspond to the descriptions presented in 
paragraphs 3.5.1 and 3*5.2. 
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FIGURE 1 

Control System Block Diagram 
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Figure 3 

Flow Chart of FINRES 
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Figure 3 (continued) 
Flow Chart of FINRES 
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Figure 3 (concluded) 
Flow Chart of FINRES 
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Figure 4 

Sample Problem Input Data 


SAHPLE PROBLEH SCOUT S-lSaC FIRST STAGE FIN RESPONSE ANALVSIS FIRST IS SEC 


1 

1 5 

2 








0. 

0. 

0. 

0 . 

0. 

0. 

0. 



0. 

0. 

0.026 0. 

0283 

0.0383 

0. 

0 . 

1 

0.0376 0 

.0376 

0.0477 0. 

0375 

0.04 0 

.04 

0.016 

I 

7 

.15 

3.86 

7.00 

3.03 

2.98 1 

.17 

15.46 

15. 


90. 

16.77 








0 . 

1. 

-1. 







1. 

0 . 

0. 







0. 

1. 

1. 







1. 

0. 

0. 








FIN 1 NON-LINEAR 

CALIBRATION ADJUSTNENTS 



- 

30. 

0. 

-10. 

0. 

10. 

0. 

30. 

1 

4 


FIN 3 NON-LINEAR 

CALIBRATION ADJUSTNENTS 



- 

30. 

0. 

-10. 

0. 

10. 

0. 

30 . 

( 

4 


FIN 3 NON-LINEAR 

CALIBRATION ADJUSTNENTS 



- 

30. 

0. 

-10. 

0. 

10. 

0. 

30 . 

( 

4 


FIN 4 NON-LINEAR 

CALIBRATION ADJUSTNENTS 




30. 

0. 

-10. 

0. 

10. 

0. 

30 . 

i 


1. 

4. 

3. PLOT SCALE FACTORS FOR 

TINE. 

FIN, AND DELTA 

-.445 

.007 

.017 

-.008 -.031 

-.117 

.076 -.883 

.028 

-.341 

.026 

-.395 

-.005 

.001 

-.008 -.006 

-.064 

.100 -.917 

.038 

-.394 

.026 

-.348 

.007 

-.018 

-.005 .016 

-.082 

.107 -.899 

.045 

-.502 

.044 

-.397 

.013 

-.021 

-.003 .014 

-.097 

.080 -.934 

.045 

-.448 

.036 . 

-.346 

.039 

-.011 

-.008 .036 

-.111 

.046 -1.021 

.038 

-.573 

.044 

-.195 

.033 

-.011 

-.008 0.000 

-.069 

.096 -1.056 

.045 

-.591 

.044 

-.145 

.035 

-.018 

-.003 0.000 

-.114 

.097 -1.073 

.045 

-.555 

.026 

-.098 

.016 

-.018 

-.010 .035 

-.107 

.074 -1.073 

.081 

-.573 

.080 

-.047 

-.031 

-.031 

-.033 .034 

-.110 

.069 -1.126 

.081 

-.627 

.116 

.004 

-.049 

-.011 

.028 .006 

-.134 

.084 -1.143 

.098 

-.627 

.098 

.055 

-.058 

-.008 

-.069 -.007 

-.147 

.095 -1.178 

.081 

-.609 

.170 

.105 

-.107 

-.008 

.272 .002 

-.108 

.147 -1.247 

.010 

-.448 

.243 

.153 

-.734 

-.034 

.604 -.041 

-.130 

.130 -1.491 

-.035 

-.072 

-.011 

.303 

-.336 

.336 

.335 -.105 

-.083 

.140 -1.735 

-.007 

.125 

-.192 

.354 

-.255 

.033 

.315 -.021 

-.121 

.181 -1.665 

-.253 

.340 

-.609 

.305 

.314 

-.099 

.528 -.016 

-.137 

.233 -1.474 

-.745 

.519 

-.773 

.355 

.355 

-.121 

.528 -.021 

-.133 

.252 -1.613 

-.833 

.483 

-.554 

.403 

-.348 

-.127 

.652 -.028 

-.086 

.315 -2.048 

-.376 

.340 

-.264 

.453 

-.342 

.145 

.776 -.087 

-.131 

.320 -2.448 

.045 

.447 

-.174 

.504 

-.493 

.229 

.827 -.133 

-.061 

.392 -3.622 

-.060 

.769 

-.391 

.555 

-.332 

.232 

.891 -.092 

-.091 

.411 -2.587 

-.534 

1.342 

-.862 

.605 

.305 

.217 

.784 -.107 

-.055 

.443 -3.500 -1.063 

1.682 

-1.134 

.653 

.261 

.004 

.685 -.082 

-.096 

.461 -2.431 -1.237 

1.915 

-1.062 

.703 

.301 

-.055 

.639 -.035 

-.101 

.483 -3.396 

-.886 

1.825 

-.637 

.754 

.091 

.014 

.558 -.045 

-.119 

.513 -3.570 

-.341 

1.664 

-.346 

.805 

-.265 

.111 

.556 -.071 

-.057 

.532 -3.796 

.116 

1.531 

-.101 

.855 

-.234 

.264 

.492 -.071 

-.109 

.546 -3.883 

.116 

1.610 

-.210 

.903 

-.033 

.339 

.436 -.065 

-.073 

.526 -2.831 

-.388 

1.915 

-.464 

.953 

.205 

.201 

.388 -.009 

-.094 

.584 -2.605 

-.640 

2.237 

-.627 

1.004 

.431 

.086 

.359 .005 

-.073 

.566 -2.361 

-.693 

2.291 

-.609 

1.055 

.357 

-.018 

.195 .100 

-.034 

.583 -2.274 

-.483 

3.165 

-.282 

1.105 

.103 

-.030 

.056 .190 

-.057 

.578 -2.336 

.045 

1.789 

.134 

1.153 

-.049 

.098 

.076 .280 

-.059 

.615 -2.431 

.625 

1.575 

.551 

1.303 

-.101 

.223 

0.000 .299 

-.006 

.610 -2.500 

1.012 

1.539 

.841 

1.354 

.044 

.261 

-.036 .446 

-.048 

.587 -2.413 

1.258 

1.736 

1.203 

1.305 

.326 

.242 

-.102 .570 

T.037 

.603 -2.135 

1.662 

1.915 

1.801 

1.355 

.404 

.104 

-.104 .711 

-.041 

.582 -1.839 

2.400 

3.023 

2.852 

1.403 

.292 

.001 

-.152 .803 

-.042 

.599 -1.717 

3.490 

1.951 

4.030 

1.453 

.035 

.017 

-.163 .872 

-.078 

.571 -1.753 

4.614 

1.700 

5.078 

1.504 

-.361 

.095 

-.127 1.030 

-.030 

.533 -1.856 

5.580 

1.557 

5.705 

1.555 

-.369 

.333 

-.096 1.074 

-.070 

.515 -1.909 

6.323 

1.575 

6.036 

1.605 

-.311 

.379 

-.081 1.129 

.036 

.562 -1.874 

6.476 

1.773 

6.351 

1.653 

-.334 

.320 

-.162 1.250 

.037 

.561 -1.717 

6.691 

1.968 

6.653 

1.703 

-.300 

.136 

-.137 1.323 

-.021 

.523 -1.491 

7.139 

2.094 

7.184 

1.754 

-.360 

.042 

-.132 1.40V 

-.062 

.493 -1.369 

7.919 

2.004 

7.875 

1.805 

-.690 

.004 

-.183 1-447 

-.024 

.563 -1.369 

8.719 

1.843 

8.437 

1.855 

-.959 

.098 

-.183 1.498 

-.019 

.514 -1.456 

9.304 

1.713 

8.743 

1.903 

-1.091 

.317 

-.239 1.537 

-.003 

.507 -1.526 

9.479 

1.682 

8.823 

1.953 

-1 . 104 

.395 

-.381 1.603 

.016 

.469 -1.439 

9.382 

1.754 

8.823 

3.004 

-1.036 

.282 

-.454 1.S48 

.043 

.465 -1.126 

9.323 

1.843 

8.855 
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Figure 4 (continued) 
Sample Problem Input Data 


s.ess - 1.079 

.179 

-.548 

1.650 

.047 

.458 

-.778 

9.362 

1.807 

9.048 

2.105 - 1 . 2 SS 

.057 

-.597 

1.704 

.005 

.366 

-.517 

9.655 

1.664 

9.305 

e.isa - 1.573 

.048 

-.597 

1 .709 

.022 

.393 

-.395 

10.084 

1.467 

9.514 

2.203 - 1.832 

. 104 

-.650 

1.681 

.059 

.389 

-.377 

10.283 

1.235 

9.546 

2.254 - 1.977 

.257 

-.645 

1.686 

.059 

.331 

-.325 

10.124 

1.145 

9.353 

2.305 - 1.968 

.348 

-.617 

1.725 

.083 

.290 

-.186 

9.733 

1.163 

9.048 

2.355 - 1.968 

.304 

-.591 

1.721 

.082 

.264 

.110 

9.343 

1.360 

8.823 

2.402 - 2.929 

.198 

-.680 

1.760 

.070 

.222 

.440 

9.109 

1.449 

8.694 

2.453 - 2.190 

.073 

-.698 

1.756 

.050 

.170 

.667 

9.070 

1.360 

8.614 

2.504 - 2.443 

.029 

-.645 

1.719 

.059 

.193 

.736 

9.031 

1.091 

8.469 

2.555 - 2.653 

.104 

-.467 

1.667 

.105 

.184 

.667 

3.836 

.948 

8.116 

2.605 - 2.752 

.195 

-.373 

1.597 

.149 

.189 

.562 

8.368 

.930 

7.585 

2.652 - 2.752 

.242 

-.228 

1.600 

.127 

.170 

.527 

7.705 

1.109 

7.055 

2.703 - 2.724 

.198 

-.051 

1.592 

.076 

.147 

.527 

7.061 

1.414 

6.653 

2.754 - 2.755 

.070 

.091 

1.611 

.062 

.141 

.545 

6.613 

1.664 

6.348 

2.805 - 2.893 

-.033 

.142 

1.528 

.120 

.197 

.493 

6.340 

1.825 

6.123 

2.855 - 3.079 

-.036 

.117 

1.443 

.111 

.202 

.249 

6.087 

1.736 

5.785 

2.902 - 3.205 

.023 

.119 

1.440 

.120 

.238 

-.012 

5.677 

1.539 

5.239 

2.953 - 3.251 

.129 

.112 

1.354 

.115 

.177 

-.151 

5.036 

1.557 

4.609 

3.004 - 3.214 

.157 

.048 

1.340 

.118 

.188 

-.134 

4.386 

1.628 

3.939 

3.055 - 3.137 

.095 

-.013 

1.260 

.135 

.208 

-.029 

3.771 

1.843 

3.432 

3.105 - 3.137 

-.005 

-.058 

1.226 

.143 

.209 

.110 

3.279 

1.843 

3.088 

3.152 - 3.208 

-.083 

-.084 

1.177 

.117 

.213 

.162 

2.963 

1.700 

2.743 

3.203 - 3.304 

-.090 

-.043 

1.111 

.143 

.218 

.127 

2.646 

1.539 

2.345 

3.254 - 3.372 

-.033 

-.051 

1.043 

.125 

.188 

.040 

2.225 

1.414 

1.783 

3.305 - 3.350 

.020 

-.063 

.967 

.115 

.217 

-.064 

1.662 

1.414 

1.203 

3.355 - 3.254 

.036 

-.109 

.956 

.125 

.182 

-.099 

.994 

1.431 

.642 

3.402 - 3.159 

-.011 

-.165 

.889 

.120 

.176 

-.064 

.485 

1.431 

.315 

3.453 - 3.094 

-.108 

-.203 

.862 

.120 

.171 

.040 

.151 

1.342 

.062 

3.504 - 3.103 

-.168 

-.226 

.815 

.087 

.171 

.075 

-.042 

1.127 

-.156 

3.555 - 3.143 

-.158 

-.251 

.815 

.053 

. 103 

.040 

-.218 

.948 

-.373 

3.605 - 3.153 

-.105 

-.284 

.722 

.099 

.123 

.005 

-.499 

.715 

-.591 

3.652 - 3.097 

-.049 

-.305 

.701 

.033 

.095 

-.012 

-.798 

.662 

- 1.007 

3.703 - 2.983 

-.040 

-.272 

.669 

.040 

.094 

-.012 

- 1.202 

.644 

- 1.315 

3.754 - 2.856 

-.105 

-.198 

.639 

.037 

.092 

.005 

- 1.554 

.662 

- 1.551 

3.805 - 2.798 

-.183 

-.129 

.566 

.081 

.126 

.023 

- 1.729 

.644 

- 1.659 

3.855 - 2.770 

-.221 

-.086 

.555 

.049 

.124 

-.047 

- 1.817 

.536 

- 1.696 

3.902 - 2.761 

-.186 

-.043 

.551 

.019 

.080 

-.238 

- 1.888 

.357 

- 1.786 

3.953 - 2.736 

-.133 

.056 

.522 

.023 

.093 

— . 464 

- 1.940 

.286 

- 1.877 

4.004 - 2.628 

-.099 

.145 

.472 

.059 

. 182 

-.586 

- 2.081 

.322 

- 2.022 

4.055 - 2.502 

-.099 

.211 

.456 

.055 

.173 

-.691 

- 2.221 

.465 

- 2.185 

4.105 - 2.394 

-.152 

.297 

.470 

-.020 

.185 

-.795 

- 2.345 

.554 

- 2.149 

4.152 - 2.335 

-.205 

.335 

.486 

-.032 

.187 

-.934 

- 2.327 

.680 

- 2.185 

4.203 - 2.295 

-.221 

.371 

.440 

.009 

.220 - 1.143 

- 2.204 

.572 

- 2.185 

4.254 - 2.267 

-.186 

.338 

.470 

.002 

.222 - 1.404 

- 2.083 

.554 

- 2.040 

4.305 - 2.215 

-.124 

.454 

.444 

-.023 

.304 - 1.595 

- 2.011 

.554 

- 1.967 

4.355 - 2.097 

-.068 

.520 

.442 

-.029 

.277 - 1.752 

- 1.975 

.662 

- 1.931 

4.402 - 1.992 

-.071 

.553 

.443 

-.001 

.348 - 1.891 

- 1.905 

.841 

- 1.804 

4.453 - 1.903 

-.127 

.596 

.475 

-.041 

.355 - 1.978 

- 1.747 

1.020 

- 1.678 

4.504 - 1.844 

-.174 

.609 

.492 

-.061 

.372 - 2.100 

- 1.431 

1.038 

- 1.424 

4.555 - 1.832 

-.193 

.650 

.510 

-.060 

.401 - 2.309 

- 1.114 

1.073 

- 1.170 

4.605 - 1.798 

-.127 

.652 

.491 

-.018 

.473 - 2.587 

-.868 

1.056 

-.935 

4.652 - 1.752 

-.033 

.627 

.502 

-.047 

.492 - 2.813 

-.657 

1.091 

-.699 

4.703 - 1.678 

.017 

.609 

.552 

-.076 

.517 - 2.900 

-.499 

1.288 

-.536 

4.754 - 1.607 

.004 

.543 

.622 

-.079 

.511 - 2.900 

-.341 

1.485 

-.319 

4.805 - 1.554 

-.049 

.447 

.615 

-.082 

.554 - 2.900 

-.042 

1.539 

-.065 

4.855 - 1.530 

-.086 

.371 

.694 

-.083 

.579 - 2.900 

.274 

1.575 

.243 

4.902 - 1.520 

-.071 

.325 

.691 

-.077 

.583 - 2.900 

.643 

1.431 

.587 

4.953 - 1.520 

-.021 

.340 

.685 

-.086 

.577 - 2.918 

.977 

1.414 

.841 

5.004 - 1.499 

.054 

.330 

.702 

-.041 

.619 - 2.970 

1.223 

1.414 

1.131 

5.055 - 1.462 

.120 

.241 

.755 

-.032 

.621 - 3.005 

1.469 

1.521 

1.348 

5.105 - 1.425 

.095 

.137 

.761 

-.079 

.618 - 2.953 

1.750 

1.557 

1.656 

5.152 - 1.406 

.048 

-.028 

.854 

-.085 

.594 - 2.796 

2.066 

1.610 

1.964 

5.203 - 1.400 

.010 

-.076 

.833 

-.058 

.613 - 2.622 

2.333 

1.521 

2.308 

5.254 - 1.406 

.007 

-.066 

.864 

-.079 

.597 - 2.535 

2^699 

1.431 

2.616 

5.305 - 1.419 

.061 

-.091 

.956 

-.093 

.591 - 2.513 

3.033 

1.396 

2.852 

5.355 - 1.437 

.129 

-.033 

.852 

-.093 

.591 - 2.552 

3.279 

1.306 

3.142 

5.402 - 1.428 

.157 

-.028 

.860 

-.093 

.590 - 2.518 

3.490 

1.396 

3.396 

5.453 - 1.422 

.161 

-.005 

.830 

-.086 

.588 - 2.448 

3.806 

1.610 

3.468 

5.504 - 1.437 

.101 

.053 

.902 

-.080 

.584 - 2.378 

3.982 

1.557 

3.993 
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Figure 4 (continued) 
Sample Problem Input Data 


S . 5 S 5 - 1.449 

.082 

.018 

1.991 

-.977 

.582 

- 2-170 

4.952 

1.664 

4.102 

5 . 6 d 5 - 1.493 

.126 

.929 

1.106 

-.953 

.591 

- 2.294 

4.386 

1.718 

4.374 

5 . 6 Sa - 1.545 

.176 

.938 

1.138 

-.074 

.622 

- 2.274 

4.632 

1.754 

4.628 

5.703 - 1.569 

.217 

.098 

1.172 

-.963 

.609 

- 2.239 

4.720 

1.843 

4.736 

5.754 - 1.576 

.211 

-.023 

1.225 

-.062 

.564 

- 2.152 

4.843 

1.933 

4.899 

5.895 - 1.585 

.161 

.010 

1.227 

-.048 

.615 

- 2.039 

4.984 

2.976 

5.944 

5.855 - 1.585 

.114 

-.098 

1.253 

-.034 

.617 

- 1.943 

5.170 

2.147 

5.297 

S . 90 a - 1.659 

.123 

.013 

1.279 

-.937 

.611 

- 1.943 

5.326 

2.139 

5.351 

5.953 - 1.737 

.173 

.051 

1.285 

-.019 

.611 

- 1.943 

5.443 

2.130 

5.464 

6.094 - 1.761 

.251 

.043 

1.268 

.018 

.579 

- 1.926 

5.560 

2.183 

5.560 

6.055 - 1.795 

.276 

.010 

1.252 

-.015 

.591 

- 1.874 

5.619 

2.326 

5.625 

6.105 - 1.789 

.239 

-.003 

1.225 

-.004 

.584 

- 1.752 

5.658 

2.452 

5.657 

6 . 1 Sa - 1.773 

.179 

-.991 

1.163 

.053 

.608 

- 1.595 

5.677 

2.488 

5.625 

6.203 - 1.832 

.142 

-.157 

1.194 

.061 

.649 

- 1.508 

5.560 

2.470 

5.480 

6.254 - 1.993 

.148 

-.157 

1.049 

.062 

.634 

- 1.474 

5.346 

2.416 

5.239 

6.305 - 1.968 

.192 

-.180 

.994 

.056 

.592 

- 1.421 

5.919 

2.416 

4.845 

6.355 - 1.989 

.254 

-.190 

9.959 

.990 

.459 

- 1.352 

4.632 

2.344 

4.356 

6.402 - 1.968 

.311 

-.058 

.931 

.072 

.538 

- 1.213 

4.122 

2.389 

3.848 

6.453 - 1.915 

.257 

-.096 

.886 

.996 

.557 

- 1 .056 

3.683 

2.452 

3.414 

6.594 - 1.847 

.239 

-.312 

.748 

.298 

-.113 

-.830 

3.349 

2.559 

3.178 

6.555 - 1.854 

.179 

-.241 

.807 

.098 

.536 

-.743 

2.910 

2.720 

2.635 

6.605 - 1.884 

.129 

-.254 

.744 

.074 

.519 

-.569 

2.629 

2.434 

2.526 

6.652 - 1.909 

.148 

-.234 

.697 

.093 

.470 

-.499 

2.418 

2.505 

2.163 

6.793 - 1.884 

.198 

-.168 

.693 

.096 

.497 

-.586 

2.066 

2.434 

1.891 

6.754 - 1.823 

.242 

-.137 

.639 

.125 

.507 

-.604 

1.680 

2.541 

1.457 

6.895 - 1.749 

.223 

-.102 

.566 

.124 

.519 

-.517 

1.311 

2.684 

1.149 

6.855 - 1.693 

.157 

-.081 

.531 

.128 

.495 

-.412 

.994 

2.792 

.950 

6.992 - 1.638 

.979 

-.084 

.524 

.158 

.552 

-.360 

.836 

2.846 

.732 

6.953 - 1.622 

.070 

-.109 

.499 

.141 

.507 

-.377 

.696 

2.828 

.623 

7.094 - 1.567 

.191 

-.119 

.443 

.171 

.493 

-.395 

.538 

2.756 

.469 

7.055 - 1.483 

.164 

-.140 

.423 

.165 

.483 

-.377 

.362 

2.667 

.315 

7.105 - 1.415 

.198 

-.203 

.412 

.175 

.473 

-.343 

.294 

2.684 

.225 

7.152 - 1.344 

.154 

-.267 

.411 

.179 

.434 

-.256 

.116 

2.720 

.207 

7.293 - 1.289 

.061 

-.294 

.377 

.173 

.448 

-.064 

.098 

2.720 

. 134 

7.254 - 1.261 

-.002 

-.330 

.399 

.156 

.394 

.040 

.098 

2 . 5:>9 

.134 

7.305 - 1.199 

-.008 

-.284 

.393 

.132 

.383 

.023 

.081 

2.416 

.134 

7.355 - 1.113 

.032 

-.211 

.359 

.169 

.390 

.023 

.063 

2.326 

.134 

7.402 - 1.054 

.070 

-.170 

.326 

.216 

.433 

.005 

.081 

2.326 

.116 

7.453 -.983 

.067 

-.190 

.358 

.192 

.362 

-.064 

.098 

2.291 

.134 

7.594 -.922 

.032 

-.162 

.381 

.145 

.389 

-.047 

.133 

2.344 

.170 

7.555 -.872 

-.021 

-.099 

.404 

.132 

.349 

-.012 

.256 

2.326 

.225 

7.605 -.817 

-.046 

-.066 

.378 

.162 

.377 

-.029 

.379 

2.326 

.334 

7.652 -.761 

-.030 

-.069 

.363 

.137 

.383 

-.134 

.520 

2.255 

.497 

7.703 -.709 

-.011 

-.020 

.421 

.149 

.385 

-.203 

.696 

2.273 

.642 

7.754 -.647 

.010 

.015 

.429 

.135 

.359 

-.238 

.871 

2.255 

.859 

7.805 -.592 

.097 

.051 

.438 

.149 

.377 

-.308 

1.047 

2.291 

1.058 

7.855 -.564 

-.005 

.074 

.463 

.157 

.407 

-.325 

1.293 

2.291 

1.258 

7.902 -.521 

-.011 

.122 

.492 

.098 

.372 

-. 325 , 

1.557 

2.326 

1.493 

7.953 -.496 

-.036 

.150 

.498 

.144 

.393 

-.430 

1.785 

2.398 

1.729 

8.094 -.471 

-.024 

.178 

.472 

.167 

.457 

-.534 

2.049 

2.389 

1.946 

8.055 -.443 

.004 

.200 

.524 

.149 

.478 

-.621 

2.277 

2.416 

2.200 

8.105 -.403 

.032 

.244 

.562 

.158 

.401 

-.708 

2.471 

2.416 

2.417 

8.152 -.379 

.029 

.294 

.697 

.111 

.404 

-.743 

2.699 

2.595 

2.671 

8.203 -.335 

.914 

.221 

.619 

.157 

.491 

-.795 

2.980 

2.702 

2.870 

8.254 -.311 

-.005 

.175 

.612 

.145 

.499 

-.812 

3.244 

2.738 

3.196 

8.395 -.329 

-.011 

-170 

.653 

.144 

.501 

-.830 

3.525 

2.720 

3.468 

8.355 -.351 

.023 

.188 

.659 

.142 

.479 

-.839 

3.806 

2.684 

3.704 

8.492 -.382 

.067 

-195 

.680 

.130 

.521 

-.917 

3.999 

2.738 

3.975 

8.453 -.373 

.104 

.170 

.740 

.115 

.498 

-.969 

4.158 

2.792 

4.120 

8.504 -.317 

.079 

.117 

.771 

.136 

.532 

-.969 

4.281 

2.899 

4.301 

8.555 -.311 

.048 

.058 

.827 

.142 

.538 

-.899 

4.474 

2.953 

4.573 

8.605 -.345 

.061 

.170 

.804 

.140 

.584 

-.830 

4.720 

2.953 

4.791 

8.652 -.354 

.073 

.183 

.828 

.159 

.545 

-.812 

4.966 

2.971 

4.999 

8.703 -.391 

.092 

.099 

.863 

.150 

.559 

-.812 

5.190 

3.042 

5.223 

8.754 -.406 

.098 

.104 

.906 

.141 

.553 

-.830 

5.326 

3.096 

5.319 

8.805 -.438 

.073 

.071 

.929 

.133 

.554 

-.847 

5.443 

3.132 

5.448 

8.855 -.453 

.945 

.048 

.936 

.139 

.524 

-.812 

5.560 

3.996 

5.641 

8.902 -.493 

.029 

.094 

.952 

.133 

.584 

-.760 

5.716 

3.150 

5.705 

8.953 -.539 

.073 

.102 

.974 

.150 

.554 

-.760 

5.833 

3.114 

5.834 

9.004 -.564 

. 101 

.063 

-961 

.174 

.598 

-.760 

5.950 

3 . 168 

5.914 
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Figure 4 (continued) 
Sample Problem Input Data 


9.055 

-.598 

.061 

.033 

1.010 

.166 

.615 

-.778 

6.088 

3.150 

5.978 

9.105 

-.686 

.048 

-.030 

1.018 

.148 

.540 

-.760 

6.126 

3.857 

6.043 

9 .isa 

-.638 

.039 

-.046 

.998 

.160 

.559 

-.760 

6.186 

3.168 

6.075 

9 . a 03 

-.650 

.089 

-.091 

.983 

.803 

.591 

-.760 

6.164 

3.150 

6.139 

9 . 85-4 

-.666 

.089 

-.190 

.999 

.173 

.555 

-.743 

6.803 

3.025 

6.187 

9.305 

-.697 

.080 

-.888 

1.069 

.158 

.516 

-.673 

6.848 

2.953 

6.219 

9.355 

-.737 

.086 

-.808 

1.088 

.130 

.496 

-.604 

6.881 

8.810 

6.835 

9.408 

-.798 

.045 

-.180 

1.063 

.181 

.509 

-.551 

6.359 

8.756 

6.284 

9.453 

-.848 

.048 

-.160 

1.068 

.178 

.541 

-.534 

6.398 

2.720 

6.300 

9.504 

-.860 

.070 

-.094 

1.056 

.164 

.589 

-.534 

6.398 

2.774 

6.300 

9.555 

-.878 

.061 

-.099 

1.055 

.175 

.480 

-.517 

6.379 

8.810 

6.300 

9.605 

-.878 

.051 

-.066 

1.079 

.198 

.534 

-.499 

6.340 

2.363 

6.884 

9.658 

-.918 

.051 

.010 

1.077 

.188 

.505 

-.464 

6.301 

2.881 

6.268 

9.703 

-.971 

.064 

.036 

1.068 

.157 

.519 

-.464 

6.301 

8.935 

6.219 

9.754 - 1.036 

.095 

.091 

1.051 

.160 

.517 

-.499 

6.831 

3.085 

6.219 

9.805 - 1.070 

.107 

.066 

1.034 

.804 

.548 

-.488 

6.848 

3.114 

6.155 

9.855 - 1.048 

.095 

.091 

1.000 

.803 

.575 

-.517 

6.186 

3.839 

5.994 

9.908 - 1.039 

.098 

.187 

1.046 

.189 

.575 

-.499 

5.970 

3.311 

5.914 

9.953 - 1.073 

.088 

. 190 

1.041 

.186 

.538 

-.499 

5.814 

3.365 

5.818 

10.004 - 1.141 

.089 

.816 

.963 

.836 

.649 

-.517 

5.736 

3.418 

5.769 

10.055 - 1.199 

.117 

.883 

.971 

.838 

.613 

-.517 

5.716 

3.490 

5.705 

10.105 - 1.818 

.145 

.341 

.978 

.839 

.687 

-.551 

5.580 

3.579 

5.496 

10.158 - 1.196 

.145 

.193 

.970 

.813 

.648 

-.551 

5.404 

3.705 

5.319 

10.803 - 1.184 

.180 

.173 

.946 

.881 

.638 

-.551 

5.190 

3.830 

5.807 

10.854 - 1.178 

.086 

. 190 

.967 

.195 

.597 

-.569 

5.001 

3.884 

5.078 

10.305 - 1.190 

.114 

.188 

.954 

.805 

.654 

-.569 

4.913 

3.908 

4.990 

10.355 - 1.818 

. 167 

.816 

.988 

.886 

.680 

-.569 

4.885 

3.955 

4.917 

10.408 - 1.833 

.811 

.811 

.877 

.834 

.668 

-.534 

4.737 

4.045 

4.887 

10.453 - 1.840 

.804 

. ISO 

.883 

.856 

.689 

-.517 

4.667 

4.842 

4.664 

10.504 - 1.809 

.145 

.135 

.873 

.861 

.708 

-.464 

4.544 

4.385 

4.519 

10.555 - 1.156 

.095 

.104 

.895 

.831 

.634 

-.343 

4.351 

4.421 

4.392 

10.605 - 1.156 

.104 

.063 

.874 

.868 

.704 

-.303 

4.845 

4.439 

4.301 

10.658 - 1.168 

.157 

.188 

.859 

.851 

.703 

-.385 

4.888 

4.421 

4.301 

10.703 - 1.159 

.170 

. 104 

.841 

.857 

.701 

-.385 

4.158 

4.439 

4.847 

10.754 - 1.138 

.139 

.051 

.843 

.874 

.700 

-.303 

4.105 

4.498 

4.138 

10.805 - 1.101 

.089 

.094 

.861 

.869 

.686 

-.838 

4.058 

4.510 

4.084 

10.855 - 1.064 

.045 

.094 

.850 

.863 

.658 

-.838 

3.999 

4.618 

4.048 

10.908 - 1.080 

.038 

.074 

.888 

.865 

.685 

-.856 

3.988 

4.588 

3.993 

10.953 

-.996 

.048 

.089 

.814 

.861 

.709 

-.856 

3.964 

4.564 

4.012 

11.004 

-.999 

.045 

.107 

.799 

.314 

.773 

-.873 

3.964 

4.510 

4.030 

11.055 - 1.005 

.038 

.118 

.778 

.301 

.762 

-.385 

3.964 

4.588 

4.018 

11.105 

-.987 

.001 

.109 

.881 

.855 

.751 

-.377 

3.964 

4.510 

4.018 

11.158 

-.897 

-.058 

.061 

.818 

.843 

.711 

-.430 

3.938 

4.498 

4.030 

11.803 

-.838 

-.077 

.051 

.801 

.865 

.761 

-.499 

3.938 

4.439 

4.030 

11.854 

-.883 

-.103 

.010 

.864 

.858 

.768 

-.569 

3.999 

4.878 

4.018 

11.305 

-.838 

-.065 

0.000 

.806 

.866 

.754 

-.656 

4.087 

4.188 

4.102 

11.355 

-.888 

-.018 

.013 

.809 

.858 

.694 

-.743 

4.810 

4.116 

4.847 

11.408 

-.916 

.004 

.003 

.888 

.845 

.763 

-.795 

4.898 

4.081 

4.865 

11.453 

-.878 

.010 

.043 

.815 

.889 

.737 

-.778 

4.316 

4.099 

4.883 

11.504 

-.883 

-.011 

.081 

.873 

.831 

.785 

-.760 

4.898 

4.152 

4.283 

11.555 

-.764 

-.058 

.099 

.831 

.836 

.768 

-.760 

4.898 

4.158 

4.883 

11.605 

-.784 

-.058 

.089 

.830 

.813 

.780 

-.795 

4.333 

4.188 

4.398 

11.658 

-.774 

-.049 

.096 

.865 

.807 

.785 

-.888 

4.481 

4.206 

4.537 

11.703 

-.789 

-.011 

.108 

.881 

.183 

.741 

-.958 

4.597 

4.158 

4.646 

11.754 

-.771 

.083 

.094 

.889 

.808 

.717 

- 1.081 

4.708 

4.188 

4.718 

11 .805 

-.764 

.083 

.094 

.885 

.834 

.816 

- 1.039 

4.708 

4.206 

4.791 

11.855 

-.746 

-.005 

.061 

.879 

.838 

.774 

- 1.081 

4.708 

4.278 

4.845 

11.908 

-.730 

-.040 

.036 

.890 

.813 

.815 

- 1.039 

4.755 

4.860 

4.936 

11.953 

-.758 

-.058 

-.085 

.908 

.817 

.787 

- 1.039 

4.878 

4.842 

4.972 

18.004 

-.798 

-.053 

-.094 

.908 

.801 

.793 

- 1.004 

5.001 

4.116 

5.044 

12 . OSS 

-. 8-18 

-.036 

-.122 

.914 

.135 

.769 

- 1.039 

5.054 

3 . 9 SS 

5 . 1-18 

18.105 

-.897 

-.011 

-.188 

.899 

.815 

.778 

- 1.081 

5.098 

3.848 

5.186 

18.158 

-.916 

-.036 

-.831 

.919 

.808 

,740 

- 1.081 

5.054 

3.830 

5.159 

18.803 

-.888 

-.086 

-.307 

.878 

.800 

.750 

- 1.004 

5.036 

3.687 

5.142 

18.854 

-.878 

-.136 

-.318 

.899 

.171 

.735 

-.958 

5.036 

3.418 

5.186 

18.305 

-.900 

-.181 

-.887 

.919 

.160 

,715 

-.917 

5.019 

3.389 

5.142 

18.355 

-.943 

-.086 

-.889 

.938 

.168 

.669 

-.899 

5.001 

3.186 

5.186 

18.408 

-.977 

-.046 

-.846 

.871 

.179 

.690 

-.958 

4.934 

3.114 

5.186 

18.453 

-.990 

-.043 

-.195 

.888 

.178 

.670 

-.969 

4.948 

3.114 

5.094 

18.504 

-.968 

-.046 

-.157 

.854 

.171 

.710 

-.969 

4.896 

3.138 

5.026 
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Figure 4 (concluded) 
Sample Problem Input Data 


12 . 5 S 5 

-.965 

-.071 

-.187 

.879 

.147 

.675 -.969 

4. 843 

3.138 

4.978 

18.605 

-.977 

-.099 

-.184 

.918 

.150 

.688 -.969 

4.773 

3.114 

4.936 

18.658 

-.977 

-.074 

-.118 

.903 

.180 

.648 - 1.081 

4.737 

3.078 

4.917 

18.703 

- 1.008 

-.049 

-.091 

.980 

.119 

.651 - 1.186 

4.755 

8,971 

4.863 

18.754 

- 1.033 

-.018 

-.074 

.884 

.144 

.643 - 1.143 

4.737 

8.935 

4.845 

18.805 

- 1.008 

.007 

-.041 

.860 

.141 

.684 - 1.160 

4.708 

3.007 

4.887 

18.855 

- 1.087 

.007 

-.015 

.860 

.187 

.686 - 1.178 

4.635 

3.060 

4.887 

18.908 

- 1.008 

.080 

.010 

.883 

.113 

.633 - 1.178 

4.614 

3.114 

4.754 

18.953 

-.968 

.086 

.061 

.871 

.131 

.657 - 1,178 

4.597 

3.138 

4.736 

13.004 

- 1.080 

.057 

.091 

.848 

.150 

.788 - 1.178 

4.568 

3.875 

4.754 

13.055 

- 1.087 

.114 

.091 

.858 

.173 

.703 - 1.160 

4.597 

3.347 

4.718 

13.105 

- 1.014 

.183 

.089 

.888 

.155 

.705 - 1.178 

4.597 

3.454 

4.664 

13.158 

- 1.005 

.101 

.096 

.859 

.148 

.681 - 1.143 

4.544 

3.544 

4.646 

13.803 

-.946 

.095 

.089 

.838 

.130 

.788 - 1.186 

4.474 

3.579 

4.591 

13.854 

-.919 

.079 

.084 

.838 

.168 

.788 - 1.108 

4.404 

3.687 

4.591 

13.305 

-.968 

.111 

.084 

.834 

.159 

.735 - 1.073 

4.404 

3.705 

4.573 

13.355 

-.968 

.195 

.114 

.850 

.174 

.783 - 1.091 

4.481 

3.741 

4.591 

13.408 

-.946 

.889 

.074 

.887 

.197 

.740 - 1.091 

4.456 

3.908 

4.609 

13.453 

-.946 

.880 

.079 

.839 

.185 

.709 -.958 

4.491 

4.045 

4.573 

13.504 

-.900 

.801 

.114 

.881 

.147 

.690 -.818 

4.481 

4.188 

4.591 

13.555 

-.857 

.188 

.086 

.864 

.181 

.733 -.743 

4.404 

4.331 

4.591 

13.605 

-.848 

.886 

.043 

.864 

.195 

.707 -.785 

4.481 

4.349 

4.609 

13.658 

-.889 

.876 

.010 

.895 

.815 

,701 -.638 

4.491 

4.481 

4.609 

13.703 

-.817 

.389 

-.083 

.847 

.859 

.739 -.517 

4.579 

4.588 

4.609 

13.754 

-.838 

.339 

-.030 

.848 

.831 

,761 -.308 

4.650 

4.636 

4.664 

13.805 

-.808 

.338 

-.003 

.846 

.881 

.774 -.099 

4.667 

4.797 

4.718 

13.855 

-.801 

.386 

.083 

.871 

.896 

.751 .058 

4.685 

4.976 

4.809 

13.908 

-.789 

.317 

.033 

.896 

.886 

,738 .145 

4,730 

5,155 

4.837 

13.953 

-.733 

.339 

.048 

.939 

.314 

.745 .814 

4.755 

5.890 

4.899 

14.004 

-.783 

.354 

.088 

.875 

.353 

.811 .884 

4.861 

5.484 

4.954 

14.055 

-.786 

.361 

.051 

.884 

.369 

.790 .440 

4.948 

5.585 

5.008 

14.105 

-.798 

.364 

.088 

.891 

.383 

.767 .545 

5.019 

5.693 

5.094 

14.158 

-.795 

.345 

.003 

.987 

.366 

.738 .667 

5.098 

5.760 

5.159 

14.803 

-.798 

.383 

-.003 

.919 

.398 

.809 .788 

5.118 

5.861 

5.807 

14.854 

-.795 

.317 

.080 

.943 

.381 

.786 .893 

5.131 

5.995 

5.855 

14.305 

-.805 

.348 

.018 

.948 

.418 

.760 .980 

5.190 

6.063 

5.887 

14.355 

-.838 

.357 

.036 

.968 

.486 

.765 1.038 

5.889 

6. 197 

5.319 

14.408 

-.845 

.376 

.038 

.935 

.466 

.301 1.108 

5.346 

6.348 

5.351 

14.453 

-.854 

.373 

.041 

.958 

.465 

.757 1.883 

5.365 

6.466 

5.351 

14.504 

-.888 

.357 

.038 

.954 

.480 

.774 1.363 

5.365 

6.600 

5.351 

14.555 

-.894 

.348 

.018 

.939 

.516 

.801 1.484 

5.365 

6.808 

5.368 

14.605 

-.918 

.317 

.008 

.914 

.533 

.888 1.606 

5.386 

6.869 

5.400 

14.658 

-.946 

.311 

.013 

.980 

.554 

.818 1.693 

5.346 

6.936 

5.384 

14.703 

-.956 

.301 

-.015 

.940 

.558 

.816 1.745 

5.365 

6.986 

5.384 

14.754 

-.968 

.898 

-.080 

.938 

.533 

.789 1.815 

5.346 

7.080 

5.368 

14.805 

-.974 

.873 

0.000 

.908 

.566 

.886 1.908 

5.887 

7.104 

5.319 

14.855 

-.971 

.854 

-.041 

.886 

.608 

.883 1.978 

5.848 

7.171 

5.319 

14.908 

-.993 

.817 

-.056 

.898 

.598 

.813 8.006 

5.809 

7.171 

5.319 

14.953 

-.999 

.189 

-.051 

.888 

.590 

.810 8.041 

5.151 

7.805 

5.351 

15.004 

- 1.005 

.167 

-.051 

.981 

.583 

.788 8.076 

5.118 

7.805 

5.871 
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r I yu re a 

Sample Problem Printed Output Data 
Individual Control Surfaces 


fAOC NO. t 


SNNPIC MOIIEN 

scour S'lsac 

rmsT 

STNOe FIN 

NESFONSe 

NNNLVSIS 

FIRST IS 

SEC 


S Ft 

N ONE 

X 

FIN T 

U 0 

I F t N 

T H R E 

E 

F X 


TINE 

COHNANO 

ACrUAl 

DELTA 

CONNAHO 

ACTUAL 

DELTA 

eONNANO 

ACTUAL 

DELTA 

CONNAND 

ACTUAL 

DELTA 

-.4S 

-.52 

-.81 

.29 

-.07 

.03 

-.09 

-.26 

-.40 

.15 

-.07 

.03 

-.09 

•.30 

-.82 

-1.07 

.24 

.14 

.05 

.09 

-.40 

-.57 

.16 

.14 

.04 

.10 

.00 

-1.13 

-1.25 

.12 

.03 

.03 

-.00 

-.56 

-.43 

-.14 

.83 

.19 

-.16 

.31 

-1.69 

-1.61 

-.08 

-.76 

-.76 

-.01 

.30 

.48 

-.18 

-.76 

-.61 

-.16 

.56 

-2.38 

•2.52 

.14 

-1.22 

-.90 

-.33 

1.59 

1.67 

-.08 

•1.22 

-1.03 

-.19 

.01 

-2.63 

-2.85 

.22 

-.31 

.08 

-.39 

1.65 

1.63 

.02 

-.31 

-.19 

-.11 

1.06 

-2.17 

-2.32 

.14 

.37 

-.08 

.45 

1.78 

1.84 

-.05 

.37 

.03 

.34 

1.31 

-1.52 

•1.93 

.41 

2.49 

2.25 

.24 

1.93 

2.00 

-.06 

2.49 

2.60 

-.11 

1.56 

-l.SS 

-1.87 

.32 

5.94 

6.42 

-.48 

1.46 

1.78 

-.30 

5.94 

6.22 

-.28 

1.81 

-1.24 

-1.43 

.19 

8.18 

9.12 

-.94 

1.54 

1.73 

-.19 

8.18 

8.65 

-.47 

2.06 

-.21 

-.59 

.38 

8.59 

9.60 

-1.00 

1.54 

1.57 

-.13 

8.59 

9.24 

-.64 

2.31 

.53 

.04 

.49 

7.57 

9.45 

-1.89 

1.18 

1.33 

-.15 

7.57 

8.89 

-1.32 

2.56 

.87 

.60 

.27 

6.46 

8.46 

-2.01 

.77 

.95 

-.18 

6.46 

7.72 

-1.26 

2.81 

.43 

.31 

.12 

4.37 

8.13 

-1.76 

1.50 

1.73 

-.23 

4.37 

5.85 

-1.47 

3.06 

.44 

.07 

.38 

1.84 

3.42 

-1.58 

1.67 

1.83 

-.16 

1.84 

3.17 

-1.33 

3.31 

.30 

-.08 

.38 

-.37 

1.17 

-1.55 

1.38 

1.43 

-.05 

-.37 

.80 

-1.17 

3.56 

.27 

.02 

.25 

-1.80 

-.43 

-1.37 

.69 

.75 

-.06 

-1.80 

-.56 

-1.24 

3.81 

-.02 

-.04 

.02 

-2.71 

-1.80 

-.91 

.19 

.53 

-.34 

-2.71 

-1.69 

-1.02 

4.06 

-.69 

-.77 

.09 

-3.12 

-2.30 

-.82 

.52 

.56 

-.03 

-3.12 

-2.17 

-.95 

4.31 

-1.64 

-1.71 

.08 

-2.62 

-1.98 

-.64 

.69 

.66 

.03 

-2.62 

-1.93 

-.69 

4.56 

-2.42 

-2.51 

.09 

-1.61 

-.93 

-.68 

1.16 

1.07 

.09 

-1.61 

-.99 

-.62 

4.81 

-2.66 

-2.90 

.24 

-.35 

.20 

-.55 

1.56 

1.55 

.02 

-.35 

.17 

-.52 

5.06 

-2.53 

-2.96 

.42 

1.02 

1.68 

-.66 

1.68 

1.56 

.12 

1.02 

1.58 

-.56 

5.31 

-2.12 

-2.54 

.42 

2.31 

3.22 

-.91 

1.34 

1.34 

-.00 

2.31 

3.07 

-.76 

5.56 

-2.04 

-2.20 

.16 

2.76 

4.29 

-1.53 

1.50 

1.71 

-.22 

2.76 

4.31 

-1.55 

5.81 

-1.73 

-1.97 

.24 

4.25 

5.12 

-.87 

1.81 

2.13 

-.33 

4.25 

5.17 

-.91 

6.06 

-1.36 

-1.78 

.42 

4.61 

5.65 

-1.04 

2.23 

2.44 

-.20 

4.61 

5.64 

-1.03 

6.31 

-.86 

-1.36 

.50 

3.22 

4.71 

-1.49 

2.47 

2.37 

.10 

3.22 

4.48 

-1.25 

6.56 

.92 

-.62 

1.54 

1.42 

2.70 

-1.29 

2.21 

2.52 

-.31 

1.42 

2.52 

-1.11 

6.81 

-.03 

-.44 

.41 

.19 

1.09 

-.90 

2.65 

2.78 

-.13 

.19 

.99 

-.81 

7.06 

.04 

-.35 

.39 

-.53 

.25 

-.78 

2.72 

2.68 

.04 

-.53 

.25 

-.78 

7.31 

.41 

.02 

.38 

-.46 

.67 

-.54 

2.24 

2.34 

-.10 

-.46 

.13 

-.60 

7.56 

.30 

-.03 

.34 

.08 

.35 

-.27 

2.25 

2.32 

-.07 

.08 

.31 

-.23 

7.81 

-.13 

-.32 

.19 

1.09 

1.23 

-.14 

2.14 

2.30 

-.16 

1.09 

1.21 

-.12 

8.06 

-.52 

-.68 

.16 

2.15 

2.43 

-.28 

2.46 

2.44 

.02 

2.15 

2.37 

-.22 

8.31 

-.63 

-.84 

.21 

3.36 

3.73 

-.37 

2.68 

2.70 

-.02 

3.36 

3.65 

-.29 

8.56 

-.59 

-.85 

.26 

4.28 

4.66 

-.38 

2.80 

2.96 

-.16 

4.28 

4.74 

-.45 

8.81 

-.52 

-.82 

.30 

6.12 

5.54 

-.41 

2.99 

3.11 

-.13 

5.12 

5.58 

-.45 

9.06 

-.54 

-.77 

.22 

5.47 

6.09 

-.63 

3.04 

3.22 

-.18 

5.47 

6.02 

-.55 

0.31 

-.17 

-.62 

.46 

5.48 

6.28 

-.80 

2.75 

2.83 

-.08 

5.48 

6.24 

-.76 

0.56 

-.12 

-.50 

.38 

5.46 

6.35 

-.89 

2.73 

2.86 

-.13 

5.46 

6.29 

-.82 

9.81 

-.21 

-.50 

.29 

5.07 

6.15 

-1.08 

3.05 

3.23 

-.17 

5.07 

6.04 

-.98 

10.06 

-.33 

-.64 

.21 

4.39 

5.61 

-1.22 

3.72 

3.58 

.14 

4.39 

5.65 

-1.16 

10.31 

-.25 

-.57 

.31 

3.84 

4.85 

-1.00 

3.89 

3.66 

-.06 

3.84 

4.93 

-1.09 

10.56 

-.04 

-.32 

.28 

3.30 

4.28 

-.98 

4.27 

4.43 

-.16 

3.30 

4.33 

-1 .03 

10.81 

.01 

-.24 

.25 

3.18 

4.02 

-.83 

4.34 

4.59 

-.25 

3.18 

4.05 

-.87 

11.06 

-.25 

-.36 

.12 

3.23 

3.97 

-.74 

4.39 

4.53 

-.13 

3.23 

4.01 

-.79 

11.31 

-.59 

-.72 

.13 

3.62 

4.18 

-.56 

3.96 

4.12 

-.18 

3.88 

4.20 

-.58 

11.86 

-.66 

-.79 

.13 

3.78 

4.33 

-.65 

3.92 

4.18 

-.27 

3.78 

4.37 

-.58 

11.81 

-.98 

-1.03 

.11 

4.23 

4.71 

-.48 

3.78 

4.86 

-.48 

4.23 

4.84 

-.61 

12.06 

-.87 

-1.03 

.16 

4.40 

6.08 

-.67 

3.63 

3.86 

-.24 

4.40 

5.14 

-.73 


PNCC NO. t 

SANPIC PN001.cn SCOUT S- 182 C FIRST STAGE FIN RCSPONSC NNALVStS FIRST IS SEC 


S 

FIN 

ONE 

* 

FIN 

T U 0 

t 

F I H 

T H R E 

E 

F 1 N 

o 

c 

R 

TIRE 

CONNAND 

ACTUAL 

DELTA 

CONNAND 

ACTUAL 

DELTA 

CONNAND 

ACTUAL 

delta 

CONNAND 

ACTUAL 

DELTA 

12.31 

-.77 

-.91 

.14 

4.16 

5.01 

-.84 

2.96 

3.20 

-.24 

4.16 

5.13 

-.97 

12.56 

-.84 

-.97 

.13 

3.83 

4.79 

-.96 

2.91 

3.11 

-.20 

3.83 

4.95 

-1.11 

12.81 

-1.04 

-1.17 

.13 

3.82 

4.68 

-.86 

2.74 

3.06 

-.32 

3.82 

4.82 

-.99 

13.06 

-1.06 

-1.17 

.11 

3.59 

4.59 

-1.00 

3.24 

3.45 

-.21 

3.59 

4.68 

-1.09 

13.31 

-.83 

-1.08 

.25 

3.56 

4.42 

-.85 

3.68 

3.76 

-.08 

3.56 

4.59 

-1.02 

13.56 

-.46 

-.72 

.26 

3.72 

4.42 

-.70 

3.98 

4.36 

-.39 

3.72 

4.60 

-.88 

13.81 

.30 

.01 

.28 

4.00 

4.68 

-.68 

4.75 

4.97 

-.22 

4.00 

4.78 

-.78 

14.06 

.80 

.52 

.28 

4.34 

5. SO 

-.66 

5.53 

5.67 

-.14 

4.34 

5.07 

-.73 

14.31 

1.31 

1.02 

.28 

4.54 

5.23 

-.69 

5.86 

6.19 

-.33 

4.54 

5.31 

-.77 

14.56 

1.85 

1.57 

.28 

4.55 

5.34 

-.79 

6.56 

6.87 

-.31 

4.55 

5.39 

-.84 

14.81 

2.22 

l.SS 

.27 

4.19 

5.26 

-1.07 

7.00 

7.16 

-.16 

4.19 

5.32 

-1.13 


31 


Figure 6 

Samole Problem Printed Output Data 
Pitch, Yaw and Roll Components 


PACiE NO. 3 PITCH, YAU, ROLL CONNANDS BASED ON MEASURED FINS 

SAMPLE PROBLEM SCOUT 3-192C FIRST STAGE FIN RESPONSE ANALYSIS FIRST 15 SEC 

TIME * COMMAND ^ACTUAL DELTA COMMAND ACTUAL DELTA COMMAND ACTUAL DELTA 


-.45 

-.07 

.03 

-.09 

-.20 

.14 

.04 

.10 

.06 

.03 

.11 

-.08 

.31 

-.76 

-.68 

-.08 

.56 

-1.22 

-.96 

- . 26 

.31 

-.31 

-.06 

- . 25 

1.06 

.37 

-.03 

.33 

1.31 

2.49 

2.43 

. 06 

1.56 

5.94 

6.32 

-.38 

1.31 

3.18 

3.89 

-.71 

2.06 

3.59 

9.42 

-.32 

2.31 

7.57 

9.17 

-1.60 

2.56 

6.46 

3.09 

-1.63 

2.81 

4.37 

5.99 

-1.62 

3.06 

1.34 
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APPENDIX A 


FORTRAN PROGRAM LISTING 


A complete FORTRAN source program listing is presented in the following 
pages. It starts with the main routine (FINRES) and is followed by the 
subroutines arranged in alphabetical order. There are a total of 182 cards in 
FINRES. The total program includiiag subroutines (less CALCOMP library) 
contains 675 cards. 


oooo 


PROGRAH FINRESdNPUT, OUTPUT, TAPE5-INPUT,TAPE6-0UTPUT, TAPES) 1 

THIS PROGRAM ANALYZES POST-FLIGHT MEASURED CONTROL SURFACE 2 

DEFLECTIONS AND AUTOPILOT PARAMETERS. IT PROUIDES AN ADJUST- 3 

MENT CAPABILITY AND DIGITAL FILTERING OF THE REDUCED TELEMETRY A 

DATA INCLUDING TIME SHIFTS. 5 

DIMENSION CDl(e0),EDl(20),CDa(2e),EDe(ae),CD3(Be),ED3(a0),CD4(B0), 6 

1 ED4(20),A(4,4),APS(4,4),FILTER(3,3),B(3),DATA(11 ),DL(10), 7 

2 NTIT(8),T(2000),Q(2000),R(2000),P(2000),TH(2000),PS(2000) 8 

3 ,PH(2000),D1(2000),D2(2000),D3(2000),D4(2000),DUN(2000) 9 

DATA NMl,NM2,NM3,NM4,Nn5,NH6,NM7/l0HFIN NO. 1 ,10HFIN NO. 2 , 10 

1 10HFIN NO. 3 ,10HFIN NO. 4 ,10HPITCH FINS, 11 

2 10HYAU FINS ,10HROLL FINS / 12 

C READ TITLE CARD 13 

READ(5,199) ( NTIT ( J ) , J» 1 , 8 ) 14 

199 FORMAT(8A10) 15 

• C READ INTEGER OPTIONS AND CONSTANTS 16 

READ(5,220) IOC, IOP,NPRT, MINT 17 

220 F0RMAT(4I5) 18 

• C READ TELEMETRY ADJUSTMENT FACTORS AND TIME FACTORS 19 

READ(5,230 ) ( DL ( I) , I » 1 , 10 ) , TCQ, TCR ,TCP, TCTH, TOPS, TCPH, 20 

1 TC1,TC2,TC3,TC4,AKTH,AKPS,AICPH,UC0 21 

UC0«6.2832»UC0 22 

230 FORMAT(8E10.3) 23 

C READ AUTOPILOT GAINS AND BREAK FREQUENCIES 24 

READ(5,230) CTH,CTHD,CPS,CPSD,CPH,CPHD,U1Q,U1R,U1P,UACT 25 

C READ FIN GAIN MATRIX 26 

DO 10 J-1,4 27 

10 READ(5,240) ( A ( J, I) , I « 1 , 3 ) 28 

240 FORMAT(4E10.3) 29 

C READ MEASURED FIN CALIBRATION ADJUSTMENT TABLES 30 

READ(5,250) NTl , ( CDl ( J ) , EDI ( J ) , J» 1 , NTl ) 31 

READ(5,250) NT2, (CD2( J ),ED2( J ), J-1,NT2 ) 32 

READ(5,250) NT3, ( CD3( J ) , ED3 ( J ), J- 1 , NT3 ) 33 

READ(5,250) NT4, (CD4( J ),ED4( J ), J-1,NT4 ) 34 

250 FORMAT(I5,/(8E10.3)) 35 



A-3 


C TEST FOR PLOT OPTION 36 

IF(IOP.EQ*e) GO TO 30 37 

C CALCOHP PLOT. READ IN SCALE FACTORS (UNITS PER INCH) 38 

READ(5,330) XSCALE . VSCALE, DSCALE 39 

IF(XSCALE.LE.0. ) IOP-0 -40 

IF(YSCALE.LE.0. ) IOP-0 -41 

IF(DSCALE.LE.0. ) IOP-0 43 

C TEST FOR PSEUDO-INUERSE CALCULATION OF PITCH, YAU, ROLL DEFLECTION 43 

30 IFdOC.EQ.O) GO TO 40 44 

C COMPUTE PSEUDO-INUERSE OF FIN GAIN MATRIX 45 

N-4 46 

M-3 47 

CALL PSEUDO(APS,A,N,H,NER) 48 

IF(NER.EQ.l ) GO TO 40 49 

IOC-0 50 

URITE(6,300) 51 

300 F0RMAT(//,5X,34HMATRIX IS SINGULAR, CANNOT COMPUTE , 53 

1 1X,14HPSEUD0-INUERSE ,/) 53 

40 NP-0 54 

REUIND 9 55 

Ml-1 56 

M3-1 57 

M3-1 58 

M4-1 59 

C BEGIN READING THE TELEMETRY DATA 60 

50 READ(5,104) ( DATA ( I ) , I - 1 , 1 1 ) 61 

104 F0RMAT(11F7.3) 63 

IF(EOF(5).NE.0) GO TO 70 63 

NP-NP+1 64 

C ADJUST TELEMETRY DATA FOR BIAS AND FIN CAL NON-LINEARITY 65 

DO 60 1-1,6 66 

60 DATA( I+l )-DATA( I+l )+DL( I) 67 

X-DATA(8) 68 

CALL TBLN(Y,X,CD1,ED1,NT1,M1) 69 

DATA(8)-DATA(8)+Y+DL(7) 70 



X-DATAO) 

CALL TBLN(Y.X,CDa,EDa,NTa,N2) 
DATA(9)-DATA(9)+V+DL(8) 

X>DATA(ie) 

CALL TBLN(Y,X,CD3,ED3,NT3,H3) 
DATA(10)-DATA(10)+Y+DL(9) 

X-DATA(U) 

CALL TBLN(Y,X,CD4,ED4,NT4,I14 ) 

DATAdl )»DATA(11 ) + Y + DL(10) 

C URITE ADJUSTED DATA ON TAPE 9 
URITE(9,104) (DATAd 
GO TO 50 
70 CONTINUE 

C READ TAPE 9 DATA INTO TABLES 
REUIND 9 
DO 80 J-1>NP 

80 READ(9d04) T( J ),Q( J ),R( J ),P( J ),TH( J ),PS( J ),PH( J ) 

• C INITIALIZE FILTER AND FILTER THESE PARAMETERS 

^ CALL FILFIL(UCO, FILTER, B,DTFIL) 

^ C COMPUTE PITCH ERROR SIGNAL AND RETURN IN 0 ARRAY 

• CALL ERSIG(FILTER,B,T,Q,TH,NP,NINT,CTH,CTHD,DTFIL, 

1 TCTH,TCQ>AKTH,U1Q,UACT) 

C COMPUTE VAU ERROR SIGNAL AND RETURN IN R ARRAY 

CALL ERSIG(FILTER,B,T,R,PS,NP,NINT,CPS,CPSD,DTFIL, 

1 TCPS,TCR,AKPS,U1R,UACT) 

C COMPUTE ROLL ERROR SIGNAL AND RETURN IN P ARRAY 

CALL ERSIG ( FILTER, B,T,P, PH, NP,NINT,CPH,CPHD,DTFIL, 

1 TCPH,TCP,AKPH,U1P,UACT) 

C COMPUTE FOUR CONTROL SURFACE DEFLECTIONS AND FILTER 

REUIND 9 
DO 100 J-l,NP 

100 READ(9,10S) T(J),Dl(J),Da(J),D3(J),D4(J) 

105 F0RMAT(F7.3,4ax,4F7.3) 

C FILTER AND ADJUST CONTROL SURFACES 

CALL FIN(FILTER,B,T,D1,NP,DTFIL,TC1,NINT) 


71 

7a 

73 

74 

75 

76 

77 

78 

79 

80 
81 
8a 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 



ooo 


CALL FIN(FILTER,B,T,Da,NP,DTFIL,TCa,NINT) 106 

CALL FIN(FILTER,B,T,D3,NP,DTFIL>TC3,NINT) 107 

CALL FIN(FILTER,B.T,D4,NP,DTFIL,TC4,NINT) 108 

THE PITCH, YAU, AND ROLL RECONSTRUCTED AND FILTERED 109 

DEFLECTIONS ARE AUAILABLE. NEXT CONFUTE THE COMMANDED 110 

INDIUIDUAL FINS AND STORE IN ARRAYS TH,PS,PH,DUM 111 

DO 110 j-i,NP iia 

TH(J)-A( 1,1 )»Q(J)+A(l,a)»R(J)+A(l,3)l:P(J) 113 

ps(j)-A(a,i)»Q(J)+A(a,a)*R(j)+A(a,3)*p(j) 114 

PH(J)»A(3,1 )»Q(J)+A(3,a)»R(J)+A(3,3)»P(J) 115 

110 DUM(J)*A(4,1 )»Q( J)+A(4,a)*R(J)+A(4,3)*P(J) 116 

C PREPARE AND PRINT OUTPUT OF INDIUIDUAL FINS 117 

NPAGE-1 118 

NLINE-59 119 

KPRT-NPRT 130 

DO 140 J»1,NP lai 

IF(NLIHE.LT.S9) GO TO 130 133 

URITE(6,305 ) NPAGE, (NTITd ), 1-1,8) 133 

NLINE-8 134 

NPAGE-HPAGE+1 135 

305 F0RMAT(lHl,aX,8HPAGE NO. ,13,// IX, 8A10, // , 9X, 16HX FIN ONE 136 

1 ,8X,16H« FIN T U 0,8X,19H* F I N T H R E E , 137 

a 5X,18H FIN FOUR , //, 3X, 6HTIME , 138 

3 4(34H COflMAND ACTUAL DELTA ) ,/) 139 

130 IF(KPRT.LT.NPRT) GO TO 140 130 

KPRT-0 131 

DLTl-TH(J)-DKJ) 133 

DLTa«PS(J)-D3(J) 133 

DLT3-PH(J)-D3( J) 134 

DLT4-DUH(J)-D4( J ) 135 

URITE(6,306) T(J),TH(J),Dl(J),DLTl,PS(J),D3(J),DLTa,PH( J),D3(J), 136 

1 DLT3,DUM( J),D4( J),DLT4 137 

NLINE-NLINE+1 138 

140 KPRT-KPRT+1 139 

306 F0RMAT(1X,13F8.3) 140 



ooo 


C TEST FOR PLOT 141 

IFdOP.EO.0) GO TO 160 142 

C CALCOMP PLOT OF INDIUIDUAL CONTROL SURFACES 143 

CALL CURUE(T,TH,D1,NP,NTIT,NN1,XSCALE,VSCALE,DSCALE) 144 

CALL CURUE(T,PS,D2,NP,NTIT,Nf12,XSCALE,YSCALE,DSCALE ) 145 

CALL CURUE(T,PH,D3,NP,NTIT,NM3,XSCALE,YSCALE,DSCALE) 146 

CALL CURUE(T,DUn,D4,NP,NTIT,NM4,XSCALE,YSCALE,DSCALE) 147 

C TEST FOR ERROR SIGNAL RECONSTRUCTION BASED ON MEASURED FINS 148 

160 IF(IOC.EQ.0) GO TO 180 149 

RECONSTRUCT PITCH, YAU, AND ROLL ERROR SIGNAL FROM MEASURED 150 

FIN DEFLECTIONS AND OUTPUT, PITCH, YAU, AND ROLL DEFLECTIONS 151 

BASED ON FINS 1 THRU 4 UILL BE IN ARRAYS TH,PS, AND PH 152 

NLINE«59 153 

DO 170 J»1,NP 154 

IF(NLINE.LT,59) GO TO 165 155 

URITE(6,207) NPAGE, (NTIT(I), 1-1,8) 156 

NLINE-8 157 

NPAGE-NPAGE+1 158 

KPRT-NPRT 159 

165 TH( J)-APS(1, 1 )*D1( J)+APS(1,2)«D2( J)+APS(1,3)»D3( J )+APS(l,4)»D4( J) 160 

PS( J)=APS(2, 1 )*D1( J)+APS(2,2)*D2(J)+APS(2,3)»D3( J)+APS(2,4)»D4( J) 161 

PH( J)=APS(3, 1 )*D1( J)+APS(3,2)*D2( J)+APS(3,3)»D3( J)+APS(3,4)*D4( J) 162 

C COMPUTE DIFFERENCES 163 

DLT1-Q( J)-TH( J) 164 

DLT2-R(J)-PS(J) 165 

DLT3»P(J)-PH( J) 166 

IF(KPRT.LT,NPRT) GO TO 170 167 

KPRT-0 168 

URITE(6,206) T(J ),Q( J),TH(J ),DLT1,R(J),PS( J),DLT2,P( J),PH( J),DLT3 169 

NLINE-NLINE+1 170 

170 KPRT-KPRT+1 171 

207 FORMATC lHl,aX,8HPAGE NO . , 13, 10X, 48HPITCH, YAU, ROLL COMMANDS BASED 172 

1 ON MEASURED FINS , // , 8A10, // , 8X, 1H», 7X, 173 

2 9HP I T C H ,7X, 1H»,9X,5HY A U ,9X, 1H*,7X,7HR 0 L L ,9X, 174 

3 1H»,/3X,6HTIME ,3(24H COMMAND ACTUAL DELTA ),/) 175 



IFdOP.EQ.0) GO TO 180 176 

C PLOT CALCOMP PLOT OF PITCH, VAU, AND ROLL FINS 177 

CALL CURUE(T,Q,TH,NP,NTIT,NnS,XSCALE,YSCALE,DSCALE ) 178 

CALL CURUE(T,R,PS,NP,NTIT,NM6,XSCALE,VSCALE,DSCALE) 179 

CALL CURUE(T,P,PH,NP,NTIT,NM7,XSCALE,YSCALE,DSCALE ) 180 

180 STOP 181 

END 188 


J 




ooo 


i 

I 

00 

I 


SUBROUTINE CURUE(T,CALC, ACT,NP,NTIT,NH,XS,YS,DS ) 

THIS SUBROUTINE PLOTS A CALCOMP TYPE PLOT OF THE COMHANDED 
ACTUAL, AND DIFFERENCE BETUEEN THE COHNANDED AND ACTUAL 
CONTROL SURFACE DEFLECTION. 

DIMENSION T(ae0e),cALC(a0e0),ACT(aeee),NTiT(8),x(a), Y(2) 

DATA X(1 ),X(a),Y(l),Y(2),NPP,ZL,2IC,SPACE,YMAX,DTHAX/5,0,6.2. 

1 a. 75, 2.75, 2, 0.07, 0.0, 0.04, 3*, 1.5/ 

DM--2.»YS 
DD--DS 

TF»T(NP)/XS + 1. 

NTF-TF 
TF-NTF 

CALL PLOTS(5HCAL22,0,4HPLOT) 

CALL PL0T(l.,6.,-3) 

CALL SYMBOL(2.,3.75,0.14,NTIT,0.,80) 

CALL SYMBOL ( 2., 3. 5, 0.1 4, NH,0., 10) 

CALL AXIS(0.,-2.,30HCONTROL SURFACE DEFLECTION DEG , 30,6 . , 90 . , DM, 
1 YS) 

CALL AXIS(0.,-2.,13HTIME SECONDS ,-13,TF,0. ,0. ,XS ) 

CALL DASH (T, ACT, NP,ZL, SPACE, ZK,XS,YS,1,TF,YHAX) 

CALL DASH (T, CALC, NP,ZL,ZK, SPACE, XS,YS,1,TF,YMAX) 

CALL PL0T(5.,2.95,3) 

PAI I P P Qt; P^ 

CALL SYnBOL(6!3i2.95,0.10,12HMEASURED FIN ,0.,12) 

CALL SYHBOL(6.3,2.75,0.10, 13HRECONSTRUCTED,0. , 13 ) 

CALL PL0T(5. ,2.75,3) 

CALL DASH ( X,Y,NPP,ZL,ZK, SPACE, l.,l.,l,TF, 4. ) 

CALL PLOT(0.,-4.,-3) 

CALL AXIS(0.,-2.,13HTIME SECONDS, -13, TF, 0. ,0. ,XS ) 

CALL AXIS(0.,-l.,a0HDIFFERENCE CALC-ACT, 20, 2. ,90. , DD,DS ) 

C CALCULATE AND PLOT THE DIFFERENCE (RECONSTRUCTED MINUS ACTUAL) 

DO 10 J-1,NP 

10 CALC(J)»(CALC(J)-ACT(J))/DS 

CALL DASH(T, CALC. NP.ZL. SPACE. ZK.XS.1..1.TF.DTMAX) 

CALL PL0T(16.,-2.,-3) 

CALL PLOT(0.,0.,999) 

RETURN 

END 
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39 


ooo 


I 

I 

VD 


I 


tDECK DASH 1 

SUBROUTINE DASH ( X, Y, NP, Z1 ,Z3, SPACE, XSCALE, YSCALE, LSYHB. XLIh, YLIH ) 3 

SYMBOLS, DASHED, DASHED-DOT LINES OR SOLID LINES UITH OR UITHOUT 3 

SYMBOLS BASED ON A SET OF SEQUENTIAL POINTS GIUEN IN 4 

THE INPUT 'X" ABSCISSA ARRAY AND THE 'Y" ORDINATE ARRAY 5 

DIMENSION X(1),Y(1) 6 

DO 10 I-1,NP 7 

XA-X(I )/XSCALE 8 

YA»Y( I )/YSCALE 9 



IF 

(ABS(XA) 

.GT.XL 

IM) 

GO TO 

10 



10 


IF 

(ABS(YA) 

.GT.YL 

IM) 

GO TO 

10 



11 


CALL PLOT ( 

XA,YA, 

3) 





13 


60 

TO 30 







13 

10 

CONTINUE 







14 

30 

IF 

(SPACE) 

330,31 

0,30 





15 


PLOT A BROKEN LINE 





16 

30 

K«0 








17 


PI3 

»1 .5708 







18 


Z = Z 

1 







19 


ZB» 

Z3 







30 


IF 

(Z3 .6T. 

0. ) GO TO 

40 




31 


ZB» 

Z1 







33 

40 

ZD» 

Z 







33 


LZ» 

0 







34 


SL» 

0. 







35 


NF» 

NP-1 







36 


DO 

300 J»l, 

NF 






37 


XA» 

X(J)/XSCALE 






38 


IF 

(ABS(XA) 

-XLIM 

.GT. 

0. ) 

GO 

TO 

300 

39 


XB« 

X( J + D/XSCALE 






30 


IF 

(ABS(XB) 

-XLIM 

.GT. 

0. ) 

GO 

TO 

300 

31 


YA» 

Y( J )/YSCALE 






33 


IF 

(ABS(YA) 

-YLIM 

.GT. 

0. ) 

GO 

TO 

300 

33 


YB- 

Y(J+1 )/YSCALE 






34 


IF 

(ABS(YB) 

-YLIM 

.GT, 

0. ) 

GO 

TO 

300 

35 
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DV-YB-YA 

DX-XB-XA 

IF (DX .NE. 0. ) GO TO 80 
IF (DY) 50,60,70 
50 TH--PI2 
GO TO 90 
60 TH-0. 

GO TO 90 
70 TH-Pia 
GO TO 90 

80 TH-ATAN(DY/DX) 

90 DX-XB-XA 
DY-VB-YA 

DZ»SQRT(DX*DX+DY*DY) 

C TEST TO SEE UHAT IS GOING ON 
IF (K) 100,180,220 
100 K»1 

SL-SPACE 

IF (DZ-SPACE) 110,120,150 
C SPACE IS LARGER THAN DZ 
110 SL-SL-DZ 

CALL PLOT (XB,YB,3) 

GO TO 300 

C NEXT POINT IS EXACTLY ONE SPACE 
120 K-0 

IF (LZ »NE. 0 ) GO TO 130 

ZD-ZB 

LZ-1 

GO TO 140 
130 ZD-Z 
LZ-0 
140 SL-0. 

CALL PLOT (XB,YB,3) 

GO TO 300 

C NEXT POINT MORE THAN ONE SPACE AUAY 
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43 

44 

45 
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47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 
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63 

64 
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66 

67 

68 

69 
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150 XA-XA+SPACE»COS(TH) 71 

VA-YA+SPACE»SIN(TH) 72 

IF (ABS(XA)-XLIM .GE. 0.) GO TO 300 73 

IF (ABS(YA)-YLIM .GE. 0.) GO TO 300 74 

K-0 75 

IF (LZ .NE. 0 ) GO TO 160 76 

ZD-ZB 77 

LZ-1 78 

GO TO 170 79 

160 ZD«Z 80 

LZ»0 81 

170 SL»0. 82 

CALL PLOT (XA,YA,3) 83 

GO TO 90 84 

C K-0 LINE BEING DRAUN ZD LENGTH NOT DRAUN RESUNE AS IS LINE STARTING 85 

180 IF (DZ-ZD) 190,200,210 86 

C LINE GOES AT LEAST TO NEXT POINT 87 

190 K-0 88 

ZD-ZD-DZ 89 

CALL PLOT (XB,YB,2) 90 

GO TO 300 91 

C LINE ENDS AT NEXT POINT 92 

200 K--1 93 

SL-SPACE 94 

ZD-0. 95 

CALL PLOT (XB,YB,2) 96 

GO TO 300 97 

C LINE ENDS BEFORE NEXT POINT 98 

210 K-1 99 

SL-SPACE 100 

XA-XA+ZD»COS(TH ) 101 

YA-YA + ZD»SIN(TH ) 103 

IF (ABS(XA )-XLIf1 .GE. 0.) GO TO 300 103 

IF (ABS<YA)-YLIH .GE. 0.) GO TO 300 104 

CALL PLOT (XA,YA,2) 105 
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ZD-e. 

GO TO 90 

C K-1 IS IN SPACE 
280 ZD-0. 

IF (DZ-SL) 230,240,270 
230 K-1 

SL-SL-DZ 

CALL PLOT (XB,YB,3) 

GO TO 300 
C SL-DZ 
240 K-0 

IF (LZ .NE. 0 ) GO TO 250 

ZD-ZB 

LZ-1 

GO TO 260 
250 ZD-Z 
LZ-0 

260 CALL PLOT (XB,YB,3) 

GO TO 300 

C SL IS LESS THAN DZ 
270 K-0 

IF (LZ .NE. 0 ) GO TO 280 

ZD-ZB 

LZ-1 

GO TO 290 
280 ZD-Z 
LZ-0 

290 XA-XA+SL»COS(TH) 

YA-YA+SL»SIN(TH) 

IF (ABS(XA)-XLin .GE. 0.) GO TO 300 
IF (ABS(YA)-YLIPl .GE. 0.) GO TO 300 
SL-0. 

CALL PLOT (XA,YA,3) 

GO TO 90 
300 CONTINUE 
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GO TO 370 

STRAIGHT LINE PLOT OPTION 
310 DO 330 J-I«NP 
XA-X(J)/XSCALE 
VA-Y(J)/YSCALE 
IF (ABS(XA)-XLin .GT. 0« ) ( 
IF (ABS(VA)-YLin . GT , 0. ) ( 
CALL PLOT (XA,YA,a) 

320 CONTINUE 
GO TO 370 

PLOT SYMBOLS ON LINE NO LINE 
330 NSH»IABS(LSYMB ) 

IF (LSYMB ,LT. 0 ) GO TO 3 
K --2 

GO TO 350 
340 K»-l 

350 DO 360 J*1>NP 
XA»X(J)/XSCALE 
YA-Y(J)/YSCALE 
IF (ABS(XA)-XLItl ,GT. 0. ) ' 
IF (ABS(YA)-YLin .GT. 0. ) i 


CALL SYMBOL (XA,YA, 
CONTINUE 

CALL PLOT (0.,0.,3) 

RETURN 

END 


(XA,YA,0.07,NSM,0.0,K) 


141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

IS NEGATIUE ' 151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 


oooo 


tDECK ERSIG 1 

SUBROUTINE ERSIG(FILTER,B,T,Q,TH,NP,NINT,CTH,CTHD,DT, 8 

1 DTT«,DTTHR, AD,U,UACT) 3 

THIS SUBROUTINE ADJUSTS THE ATTITUDE RATE AND DISPLACEHENT A 

FOR TELEMETRY CROSS COUPLING, TIME LAGS, AND THEN COMPUTES 5 

THE ERROR SIGNAL. THE DISPLACEMENT AND RATE ARE ENTERED 6 

IN ARRAYS TH AND Q. THE ERROR SIGNAL IS RETURNED IN (Q). 7 

DIMENSION FILTER(3,3 ),B( 3), T( 8000), 0(8000 ),TH( 8000), DUM (8000), 8 

1 A(3,3),AB(3) 9 

C FILTER THE DISPLACEMENT TIME HISTORY 10 

N=*3 11 

K-3 12 

CALL TRESP( FILTER, B,T,TH,DUM,NP,N,K, MINT) 14 

CALL TRESP(FILTER,B,T,Q,TH,NP,N,K,NINT) 15 

M»1 16 

DTT»DT+DTTM 17 

DTR-DT+DTTMR 18 

C AT THIS POINT THE FILTERED DISPLACEMENT IS IN (DUM) AND THE 19 

C FILTERED RATE IS IN (TH) 20 

DO 10 J-1,NP 21 

TA-T(J)+DTT 28 

CALL TBLN(DISP,TA,T,DUM,HP,M) 23 

TB«T(J)+DTR 24 

CALL TBLN(RATE,TB,T,TH,NP,M) 25 

C COMPUTE THE ERROR SIGNAL AND PLACE IN (Q) 26 

10 Q(J)-CTH«DISP+CT*RATE 27 

C NEXT COMPUTE THE FIRST ORDER AUTOPILOT LAG (U) AND THE FIRST 28 

C ORDER ACTUATOR LAG (UACT) RESPONSE. 29 

AB(8)-0. 30 

A(l,8)-0. 31 

N-8 38 

K-8 33 

A6(l)-U 34 

A(l,l)--U 35 

A(8,1)-UACT 36 

A(8,8)--UACT 37 

CALL TRESP(A,AB,T,Q,DUM,HP,N,K,NINT) 38 

C PUT RESULTANT ERROR SIGNAL IN (Q) ARRAY. 39 

DO 30 J-1,NP 40 

30 Q(J)-DUM(J) 41 

RETURN -<2 

END 43 



o o o 


I 





I 


»DECK FILFIL ^ ^ 

SUBROUTINE FILFILdICO, A,B,DT ) 

THIS SUBROUTINE FILLS THE FILTER COEFFICIENTS FOR A THIRD 
ORDER BUTTERUORTH FILTER HAUING A CUTOFF FREQUENCY OF 
(UCO) RADIANS PER SECOND. 

DIMENSION A(3,3),B(3) 

C ZERO OUT THE COEFFICIENT ARRAYS 

DO 10 J-1.3 
B( J )-0. 

DO 10 K«l,3 
10 A(J,K)»0. 

C FILL REMAINING CONSTANTS 

TAU«0.7071/UCO 


A(3,l )-1.0 

A(3,3)-1.0 

A(l,l )«-a./TAU 

A(1,2)«A(1.1 )/TAU 

A(1,3)--1./(TAU*TAU*TAU) 

B(1 )»-A(1.3) 

C COMPUTE THE LOU FREQUENCY TIME LAG FOR THIS FILTER 
DT»H.*TAU 
RETURN 
END 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 
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»DECK FIN 

SUBROUTINE FIN ( FILTER, B, T, D, NP, DTFIL, TC, NINT ) 

C THIS SUBROUTINE FILTERS THE CONTROL SURFACE DEFLECTION AND 
C ADJUSTS FOR TIME DELAYS. 

DIMENSION FILTER(3,3 ),B(3),T(2000),D(2000),DUn(S009) 

N-3 

K-3 


DT-DTFIL+TC 

CALL TRESP (FILTER, B,T,D,DUM,NP,N,K,NINT) 
M» 1 

C THE FILTERED DEFLECTION IS NOU IN (DUM) 
DO 10 J«1,NP 
TA»T( J ) + DT 

CALL TBLN(FN,TA,T,DUM,NP,M) 

10 D(J)»FN 
RETURN 
END 


1 
1 
12 
13 
iA 

15 

16 
17 
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XDECK 

PSEUDO 

1 



SUBROUTINE PSEUDO ( B, A, N, N. NER ) 

3 

C 


THIS SUBROUTINE COMPUTES THE PSEUDO INUERSE MATRIX B FROM THE 3 

C 


N BY M MATRIX (A). (B )-A#-INB(AT*A).AT 
DIMENSION AS(4,4),A(4,4),B(-4,4),AINU(4,4) 

4 



L-M 

5 



LS-N 

6 



IF(N.GT.M) LS-n 

7 



IF(N.GT.M) L-N 

8 

C 


SET MATRIX ELEMENTS TO ZERO 

9 



DO 10 J-1,L 

10 



DO 10 K-UL 

11 



B( J,K)»0. 

13 


10 

AS( J,K)«0. 

13 

C 


COMPUTE THE TRANSPOSE OF A 

14 



DO 30 J-1,L 

15 



DO 30 K-1,L 

16 


B0 

AS( J,K)-A(K,J ) 

17 

C 


MULTIPLY A TRANSPOSE TIMES A, AND STORE IN B 

18 



CALL XMULT(AS,A,B,L) 

19 

C 


COMPUTE INUERSE OF B AND STORE IN AINU 

30 



CALL SIMEQ(B>LS,AINU>NER) 

31 



IF(NER.NE.0) GO TO 30 

33 

C 


MATRIX IS SINGULAR RETURN TO CALLING ROUTINE 

33 



RETURN 

34 

C 


COMPUTE THE PSEUDO INUERSE 

35 


30 

CALL XMULT(AINU,AS«B«L) 

36 



RETURN 

37 



END 

38 
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>DECIC RUNGE 

SUBROUTINE RUNGE (N,F, H, X, V, L,I) 

C THIS SUBROUTINE PERFORMS THE RUNGE-KUTTA INTEGRATION 
C UPDATES FOR THE TRESP SUBROUTINE. 

DIMENSION V(1 ),F(1 ),SU(3),FF(3) 

I - I + 1 

GO TO ( 1, a. 3, 4, 5) , I 
1 L • 1 
RETURN 

a DO 6 JM,N 
SU(J) « Y(J) 

FF(J) - F(J) 

6 V(J) - SU(J) + .5*H*F(J) 

X » X + .5»H 

L = 1 
RETURN 

3 DO 7 J»1,N 

FF(J) - FF(J) + a.*F(J) 

7 Y(J) - SU(J) + ,5*H*F(J) 

L » 1 

RETURN 

4 DO 8 J*1.N 

FF(J) = FF(J) + a.»F(J) 

8 Y(J) » SU(J) + H*F(J) 

X = X + ,5*H 

L “ 1 
RETURN 

5 DO 9 J»1,N 

9 Y(J) « SU(J) + (H/6. )»(FF(J) + F(J)) 

L » a 

I » 0 

RETURN 

END 


1 

a 

3 

4 

5 

6 

7 

8 
9 

10 

11 

la 

13 

14 

15 

16 

17 

18 
19 

a0 

ai 

aa 

E3 

a4 

as 

86 

37 

as 

39 

30 

31 

38 
33 



oo 


I 


I 


I 


tDECK SINEQ (A,KC,AINU, lERR) 

DIMENSION A(4?4),B(4.4).XD0T(S),X(5),AINU(4,4) 

?Hll SUBROUTINE HNOS THE INUERSE OF MATRIX (A) USING 
DIAGONALI2ATION PROCEDURES. 


N-1 

lERR-i 

DO 10 I«1#KC 
DO 10 J-1>KC 
AINU(I,J)-0. 
10 B( I, 

DO a0 I-1,KC 
AINU(I,I)“1. 
20 X( I )»XDOT( I ) 
DO 110 I»1»KC 
COMP»0. 


30 If\aBS(B(K,I ) )-ABS(C0PIP) .LE. 0.) GO TO 40 
C0MP»B(1C,I ) 


N»K 


40 K»K+1 

IF (K-KC .LE. 0 ) GO TO 30 
IF (B(N, I ) .EG. 0. ) GO TO 120 
IF (N-I) 120,70,50 
50 DO 60 n»l,KC 
TEMP»B(I,H) 

B(I,M)»B(N,I1) 

B(N,M)-TEI1P 

TEMP»AINU(I,H) 

AINU(I,M)«AINU(N,M) 

60 AINU(N,M)»TEMP 
TEMP«X(I ) 

X(I)-X(N) 

X(N)-TEMP 

70 X(I )-X(I )/B(I,I ) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 
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21 
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24 

25 

26 
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28 
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32 

33 

34 
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TEMP-B(I,I ) 

00 80 n-l«KC 

AINU(I,H)-AINU(I,M)/TEnP 
80 B(I,«)-B(I,fl)/TEMP 
00 100 J-1,KC 

IF (J-I .EG. 0 ) GO TO 100 
IF (B(J»I) .E0« 0.) GO TO 100 
X(J)-X(J)-B(J,I)»X(I) 

TEHP-B(J,I ) 

00 90 N-l^KC 

AINV(J,N>-AINU(J,N)-TEHP»AINU(I,N) 

90 B(J,N)-B(J,N)-TEMP»B(I,N) 

100 CONTINUE 
110 CONTINUE 
RETURN 

130 URITEC 6,130) 

IERR*0 

RETURN 

130 FORMAT (6X,aaHTHE MATRIX IS SINGULAR) 
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»DECK TBLN 1 

SUBROUTINE TBLN ( Y,X,T, A>NT,H ) 8 

THIS SUBROUTINE IS A TABLE LOOKUP FROM ABSCISSA TABLE ^T^ 3 

AND ORDINATE TABLE'A^. IS ORDINATE AT GIUEN ABSCISSA 'X'. A 

^NT^ IS LENGTH OF TABLES^T'AND^A' . 'K* IS LOCATION OF LAST UALUE 5 

DIMENSION T(1),A(1) 6 

10 IF (T(M)-X) 50,20,30 7 

20 Y-A(M) 8 

RETURN 9 

30 IF (T(1)-X.LT.0. ) GO TO 40 10 

n-1 11 

GO TO 20 la 

40 n=M-l 13 

GO TO 10 14 

50 MM«M+1 15 

IF (MM-NT.LE.0) GO TO 60 16 

M = NT 17 

r GO TO 20 18 

60 IF (T(MM)-X.GT.0. ) GO TO 70 19 

n-m 20 

GO TO 50 21 

70 n»Mn-l 22 

DT»T(MM)-T(M) 23 

IF (DT.NE.0. ) GO TO 80 24 

Y«A(M) 25 

RETURN 26 

80 DY«A(MM)-A(n) 27 

DDT-X-T(M) 28 

Y«A(M)+DY»DDT/DT 29 

RETURN 30 

END 31 



ooo 


I 

I 

ro 

ro 

I 


OF 

EQUATIONS. 


SiHENSION R(l3»%<3),T(am),V(Ee00).Z(2«««).Xm(3 
STP=NINT 
Z(l)-Y(l) 


n=i 

DO 1 J-1,3 
XD0T( J)-0. 

1 X(j)»0. 

DO BO J-B,NP 
TA»T( J-1 ) 

DT=T( J)-TA 
H»DT/STP 
DO 10 I-l.NINT 
II«0 

5 CALL RUNGE(N,XDOT,H>TA.X,L,II) 
IF(L.EQ.B) GO TO 10 
CALL TBLN(U,TA,T,Y,NP,H) 

CALL VD0T(A,X,XD0T,B,U,N) 

GO TO 5 
10 CONTINUE 
Z(J)-X(K) 

B0 CONTINUE 
RETURN 
END 


3) 
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tDECIC XnULT 

SUBROUTINE XNULK A.B.C.N ) 

C THIS SUBROUTINE COMPUTES THE PRODUCT OF TUO MATRICES A AND B 

C HAUING DIMENSIONS N BY N. THE RESULT IS C 

DIMENSION A(4,4),B(4.4),C(4,4) 

C SET UORKING MATRIX ELEMENTS TO ZERO 
DO 10 J-1«N 
DO 10 K-1*N 
10 C(J,K)»0. 

DO E0 J-1,N 
DO 30 K-1,N 
DO 30 JiC-l«N 

30 C(J,K)«C(J,K)+A(J,JK)»B( JK,K) 

RETURN 

END 
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13 
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»DECK VDOT 1 

SUBROUTINE VDOT ( A, Y, XDOT» B, U, N ) 3 

C THIS SUBROUTINE UPDATES THE STATE UARIABLE EQUATIONS FOR 3 

C THE TRESP RUNGE-KUTTA INTEGRATION SUBROUTINE. 4 

DIMENSION Y(3),A(3,3),XDOT(3),B(3) 5 

DO 3 I«1,N 6 

XDOT(I)«0. 7 

DO 1 J-1,N 8 

XDOTCI )-XDOT(I )+A(I,J)*Y(J) 9 

1 CONTINUE 10 

XDOTd )-XDOT(I )4^B(I )«U 11 

3 CONTINUE 13 

RETURN 13 

END 14 
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